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ABSTRACT 



In his Ph.D. thesis Doctor Irvin R, Lindemuth, at the 
University of California/Livermore, developed a very 
general method for numerically solving two-dimensional, 
two-fluid magnetohydordynamic equations, A copy of his 
computer code was given to the Physics Department at the 
Naval Postgraduate School for conversion to the IBM 360/67 
system presently in operation at the school. This paper 
is intended to be a users manual for this code. Numerous 
changes to the original code were required due to the 
inherent differences between the CDC and IBM machines. 

The conversion of this code as well as a complete under- 
standing of its operation and logic was the goal in the 
preparation of this paper. 
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I. INTRODUCTION 



A. BACKGROUND 

In his Ph.D. thesis, THE ALTERNATING-DIRECTION IMPLICIT 
NUMERICAL SOLUTION OF TIME-DEPENDENT, TWO-DIMENSIONAL, 
TWO-FLUID MAGNETOHYDRODYNAMIC EQUATIONS, Doctor Irvin R. 
Lindemuth developed a very general and complex computer 
code which he made available to the physics department at 
the Naval Postgraduate School. The code was originally 
written in CDC FORTRAN for use on the CDC 3^00 computer. 

This machine has a relatively small magnetic core which 
makes available to the user only about 28,000 words of 
memory. Procedure modifications required to conserve core 
storage and the generality of the code contributed to its 
complexity. With the proper selection of input parameters 
the code can solve one or two dimensional problems in 
either rectagular, cylindrical, or spherical coordinates. 

It also provides extensive control of data transfer between 
the CPU and its associated 1/0 devices. Diagnostic 
subroutines are available to assist in detecting numerical 
instabilities and inaccuracies during execution of the code. 
An additional feature is the restart capability which allows 
the program to pick up execution at the point where it was 
terminated or at any previous time step which was stored 
in external memory. 

The code was obtained in the form of a deck of punched 
cards. This deck consisted of approximately 5500 source 
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cards. A test program, sample output, Lindemuth's thesis, 
and a brief description of the code were obtained at the 
same time. The author and his advisor, Doctor G. A. 

Garre ttson, consulted with Doctor Lindemuth on several 
occasions to discuss the complexities and subtilities of 
the code. With these sources of data and information, the 
conversion process was begun. 

This code is of particular interest to the Physics 
Department at the Naval Postgraduate School where projects 
researching laser produced plasmas confined in a magnetic 
field and theta-pinch plasma devices are currently under 
way. It has verified the shell structure of a laser- 
produced plasma observed by Schwirzke and should be a 
valuable tool in future research. For a detailed explana- 
tion of the equations and numerical methods incorporated 
into this code, the reader is referred to Lindemuth's 
thesis (UCRL 51103). 

B. PURPOSE 

This thesis work was strictly concerned with converting 
Lindemuth's code to a form which could be used on the 
school's IBM 360/67 computer. As a result this paper is 
intended to be a users manual for the code and will be 
devoted primarily to the construction of the code along 
with detailed descriptions of subroutines, files, arrays, 
variables, input parameters, and so forth. 
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C. ORGANIZATION OF THE PAPER 



Consistant with the subject, this paper is organized 
like a computer code which has a small main program and 
many subroutines. The main body of the thesis consists of 
short sections sketching the plan of attack for the conver- 
sion, a few comments on problems encountered, operation and 
use of the code, and a few conclusions regarding the 
conversion. The bulk of the presentation is contained in 
the appendices. These appendices contain detailed descrip- 
tions of particular elements such as arrays, subroutines, 
common blocks, input parameters, and so forth. 

II. CONVERSION PROCEDURE 



A. PLAN 

Since the code was written in FORTRAN programming 
language, it was expected to be reasonably compatable with 
the IBM machine. The plan for conversion was to place the 
source code in the IBM FORTRAN compiler and see what 
syntax differences existed. These differences would be 
resolved and the source code resubmitted to the compiler. 
This process would continue until the entire program could 
be compiled without error diagnostics. At this point the 
program would be loaded and execution would be attempted 
with the test problem input furnished by Lindemuth. 
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B. PROBLEMS ENCOUNTERED DURING COMPILATION 

After the first compilation attempt, it became apparent 
that the conversion would be considerably more complicated 
and time consuming than anticipated. The difficulties 
divided themselves into two categories. One was differences 
in the FORTRAN programming language between the CDC and IBM 
machines and the other was limitations in the IBM catalogued 
procedures for compilation. 

1 • Programming Language Differences 

The variations in programming languages were 
apparent in the source listing output of the compiler. A 
detailed list of the language differences and the action 
taken to resolve them is presented in appendix A. 

2 . Procedure Limitations 

Most of the time spent in obtaining a syntactically 
correct source deck was due to turn-around time and non- 
productive runs caused by unforseen limitations in the 
catalogued procedures used to compile the source deck. 
Appendix B discusses Job Control Language (JCL) cards and 
catalogued procedures in detail. The majority of the prob- 
lems were the direct result of the enormous size of the 
code. Lindemuth's source code contained close to 6000 cards 
after the multiple command cards were separated into single 
cards. Working with individual subroutines eliminated 
most of the size limitation problems during compilation. 
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Another size limitation was experienced during 
the compilation of subroutine MAT2 . This subroutine was 
so large it exceeded the CPU core allocated by the catalogued 
procedure. An over-ride statement allocated additional core 
and the compilation then continued to completion. 

C. PROBLEMS ENCOUNTERED DURING EXECUTION 

After the code was successfully compiled an execution 
attempt was made. At that point several additional prob- 
lems arose. 

1. Core Limitation 

At first the program failed to load into core be- 
cause the catalogued procedure for execution limits the 
program and all external references to 100K bytes of core 
memory. A check of the linkage editor output showed that 
the code required over 200K bytes of memory to load. This 
resulted in a non-productive run, but was easily corrected 
by over-riding the catalogued procedure's REGION declaration. 

2. Subroutine Calls 

It was discovered during a subsequent execution 
attempt that call were being made to subroutines which 
could not be located by the linkage editor. It was at this 
point that an additional discrepancy in programming lan- 
guages was discovered. CDC FORTRAN uses different 
identifiers for the subroutines in the standard Scientific 
Subroutine Package (SSP) which forms a portion of the soft- 
ware accompanying the machine. For example, IBM FORTRAN 
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would call EXP(arg), while CDC FORTRAN would call EXPF(arg). 
After all of the subroutine identifiers had been corrected 
the code compiled and link edited properly with a REGION 
requirement of 214 k bytes, 

3. Data Sets 

Lindemuth's code originally provided for creation 
and destruction of files, or external data sets, by calling 
a utility subroutine CREATE which is available in CDC 
FORTRAN, Appendix C discusses files and data sets in de- 
tail, The IBM system requires that the data sets be declared 
and defined before execution of the program. Therefore the 
number, size, and description of data sets had to be 
determined and entered before the program could be executed. 
All coding which involved the CREATE subroutine was removed 
from the code. 

4. Abnormal Ends (ABEND) 

During another execution attempt an "underflow" 

ABEND was received in subroutine TCINIT. This again was 
the result of machine differences. The CDC machine can 
handle floating point numbers in the range 10 + 3 22 to 10” 2 93 
while the IBM machine will only handle number to the order 
Of 10* ? 5. In this subroutine a constant was being evaluated 
which involved the multiplication of several very small 
numbers followed by division of some very small numbers. 
Although the final result was well -within the IBM limit, at 
certain points of the calculation the accumulation register 



11 



was required to store a number whose value was outside the 
limit and the ABEND resulted. The IBM control program has 
a STANDARD FIXUP subroutine which assigns a precision zero 
to the register which underflows and then continues execu- 
tion of the program. This, of course, is unacceptable in 
these calculations and frequently leads to another ABEND 
when division by this constant is attempted. This problem 
was solved by rewriting the FORTRAN statement in such a way 
that the accumulation register never exceeds the limits 
during the calculation. 

5. Typographical Errors 

This is perhaps the most elusive error encountered. 

A tremendous amount of rewriting was done to the original 
code. During these revisions a number of typographical 
errors occurred. Some of these were discovered while study- 
ing a section of the code on a subject completely unrelated 
to the error. This type of discovery is fortunate but 
largely a matter of luck. These errors have usually been 
discovered while looking for the cause of an ABEND. For 
example, during an execution attempt the program terminated 
with a completion code 0C4 which indicates that there is a 
subscripting error where the program was ABENDED. However 
this message by itself gives no indication of where in the 
program the execution terminated. The program was run 
again with a SYSUDUMP included in the JCL cards. Upon an 
ABEND the contents of core and all of the registers are 
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printed out. Among many other things the SYSUDUMP output 
contains a program source work (FSW) which gives the location 
of the command which caused the ABEND. The PSW and the 
linkage editor map can then locate the troublesome sub- 
routine. Examination of the subroutine revealed the error 
in an EQUIVALENCE statement where NADDA, an index for a 
subscripted variable, had been typed MADDA. Since MADDA 
had not been initialized in the subroutine, the computer 
picked whatever value was in its storage location and 
attempted to use it as the index. It was highly improbable 
that this value was an integer and even less probable that 
it was within the range of the array dimension. Therefore 
the program naturally ABENDED. 

III. PROGRAM ORGANIZATION 
A. ORIGINAL CODE 

The original code consists of the main program and 

thirty-three subroutines. These subroutines can be placed 

in the following categories: 

Initialization 
Logic and Control 
Operation 
Output 
Diagnostic 

These categories are discussed in appendix D which also 
contains a detailed explanation of each subroutine. 

Execution options in the program are determined by 
parameters entered on input cards. Appendix E contains 
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a list of all the input cards; this list describes the 
card format, calling routine and statement number, and 
variable identification. 

Data is passed between subroutines by the use of COMMON 
blocks of storage which are then equivalenced to operating 
variables in the subroutines. Appendix ? contains a list 
and description of the COMMON blocks used in this program. 

The computed values of the dependent variables are 
stored in the VAR vector of COMMON/C 3/. This vector to- 
gether with the DADT vector, which contains time step 
information, is buffered out to external permanent storage 
after each time step. The DA vector and the HZ vector 
which contain startup information and mesh point location, 
respectively, are buffered out initially to provide a 
restart capability at any past time step in the program ex- 
ecution. Appendix G contains a list and detailed description 
of all the important vectors and arrays used by this program. 

The CDC version of this code provides for some external 
control during program execution by using ’’sense switches.” 
These are actual switches located on the operator's console 
that flag a preset location which can be tested by a log- 
ical GO TO command. Thus some degree of control can be 
exercised during execution of the program. These switches 
are used in conjunction with subroutines CHECK and SPECHK 
to get convenient diagnostics for checking numerical 
instabilities . 
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The code is set up to buffer out data to storage after 
each time step and to print the data after each nth time 
step as determined by the input parameter NDTPNT. 

B. PRESENT CODE 

The present code consists of the main program and 
thirty-four subroutines. The organization is essentially 
the same as the original except that all sense switch 
capability has been removed. Use of the diagnostic sub- 
routines must be determined before program execution and 
requests entered on the input cards. A modification has 
been added to allow the code to buffer out data to storage 
after each mth time step as determined by the input param- 
eter NDTBUP • The additional subroutine was added to 
handle binary data transfer in and out of external storage. 
The converted code is listed beginning on page 105. 

IV. USE OF THE CODE 



A . PRESENT USE 

During the testing and initial run phases the code has 
been run using punched cards as input to OS. The source code 
was compiled to produce an object deck which was then used 
to input the program to the system 360 computer. Using 
the object deck as input has the advantages of reading in 
less cards to start the program and requiring less set-up 
time (compilation and link editing of source code) before 
execution. For example the compilation step takes 
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approximately three and one-half minutes of CPU time 
while the link step takes only four seconds. 

The test program was run using a JCL deck of the 



form 



//YAG0833 JOB ( 0833t 0739FP, FA11 ) , 'THESIS’ 
// EXEC FORTLG 

//link.sysin dd * 

(set of object decks) 

/* 

//GO.FT06F001 DD SPACE= (CYL, (3 . 1 ) ) 
//GO.FTlOFOOl DD UNIT=23l4 f DSNAME=&REC0RD, 
DCB=(RECFM=FB,LRECL=100,BLKSIZE==3500) f 
V 0L= SE R= MA R Y , SPACE= ( TRK , ( 1 5 . 1 ) . HLSE ) , 
DISP=( NEW, DELETE ) 

//GO.SYSIN DD * 

(test program input deck) 

/* 



C 

C 

c 



The JCL program format is discussed in appendix B. The 
data set definition card which begins //GO.FTlOFOOl and 
the data set over-ride card which begins //GO .FT06F001 
are discussed in detail in appendix H. 

All program runs will use this same general format. 
Program modification will be performed via the input deck. 
Appendix E lists all of the input cards and describes them 
completely. All of the input cards listed must be included 
in the input stream in the listed order, even though the 
parameters may all be zero. Two exceptions to this rule are 
also discussed in appendix E. 

To determine the input parameters the problem must be 
formulated and initial and boundary conditions fixed. The 
appropriate section of reference 1 (for example, application 
to laser produced plasmas, application to theta pinch) will 
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assist the operator in determining what input parameters 
are required. After the problem is formulated, use of 
appendix E in conjunction with appendix D (subroutine de- 
scriptions) will determine which input cards and parameters 
to include in the input data set stream. 

The output format is composed of two parts - the input 
data and the computed values of the dependent variables. 

The input data listing is a labeled listing of all the 
variables, boundary conditions, and initial conditions 
which are read in the input data set. The computed values 
are separated into sections by dependent variable. Each 
section contains a plot of the dependent variable in the 
Z plane following by a matrix listing of the value of the 
dependent variable at each mesh point. This last section 
is printed each time 0UTPT1 is called as determined by NDTPNT. 

B . FUTURE USE 

Once the code has been completely tested and is oper- 
ating satisfactorily it will be placed on a permanent 
external storage device. Methods for this storage are 
discussed in reference 8. This will have the advantages of 
less handling of card decks and increased access time and 
flexability. When the code is placed on external storage, 
it can be stored in a form which can be loaded directly into 
core without compilation or link editing. Thus by using 
only a few cards in the JCL program, the entire program code 
can be entered and made ready for execution in less than a 
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second. The only additional cards to be manually entered 
would be the input data cards which would be determined 
as before. 



V. TEST RUNS AND RESULTS 

At the time of this writing the test program had only 
partially executed successfully. The converted code did 
successfully execute the initializing functions and sub- 
routines which set up initial conditions for all of the 
variables at each mesh point and established proper 
boundary conditions. The execution ran into difficulty 
when data transfer into external storage was attempted. It 
is anticipated that this problem can be quickly solved and 
the program executed through normal completion. 

A printout of the test program output was furnished 
with the original code. A copy of this printout is included 
as part of reference 9 which has been given to Doctor 
G. A. Garretts on of the Physics Department at the Naval 
Postgraduate School. In addition to the test results, this 
printout contains notes explaining the contributions of 
specific subroutines to the output format. Reference 9 
also contains a listing of the original code, the initial 
test results of the converted program, and a compilation 
listing of the converted program which includes an object 
code listing for debugging operations. 
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VI. ERROR MESSAGES AND DIAGNOSTICS 



A complete treatment of error messages and diagnostics 
is beyond the scope of this paper. Reference 10 contains 
detailed discussions of all the error messages and completion 
codes, their meanings, probable causes for the message, and 
possible corrective actions. 

Reference 9 contains notes on the use of object code 
listings and SYSDUMP (a core dump upon an ABEND) printouts 
to locate a statement which is causing an error. To get 
a SYSDUMP insert the JCL card //GO. SYSUDUMP DD SYS0UT=A 
immediately before the //GO.SYSIN DD * card in any of 
the execution programs. The single most important piece of 
information to be extracted from the dump is the program 
source word (PSW). The last six characters of this word 
give the location address (in hexadecimal) in core of the 
command which caused the ABEND. Subtracting the entry address 
of the program from this ABEND address will give the location 
in the program of the bad command. This reference address 
can be used with the control section map in the linkage 
editor output to locate the subroutine which contains the 
bad command. To locate the exact statement which is causing 
the error subtract the origin of the subroutine from the 
reference address. This will give the address (in hexa- 
decimal) of the bad command in its subroutine. Now by 
locating this address in the object code listing, the 
statement number of the bad command can be extracted. 
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The computer center provides duty programmers to assist 
the operator in interpreting error messages which he cannot 
understand. They will explain generally what is causing the 
error and refer the operator to a specific reference which 
deals with the problem. If an error cannot readily be 
located, much time will be saved by consulting these pro- 
grammers immediately. 

VII. CONCLUSIONS 
A. PROGRAM CONVERSIONS 

Conversion of computer codes can be an extremely com- 
plex and time consuming task. Often it is a difficult task 
to convert a code which is to be used on the same type of 
machine in different installations. When the machines are 
of different manufacture the complications can reach 
enormous proportions, even though the code is written in 
the same programming language. It is frequently more 
advantageous to write a new code than to convert the old 
one . 

A good flow chart for a code is an invaluable tool. 
Without it, the programmer may spend excessive time trying 
to understand an operation which would be obvious if he 
had a flow chart which showed the logic of the operation. 
Plenty of comment cards are beneficial, especially if the 
code is to be used by others. 
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The decision to convert a computer code should be made 
carefully. The fact that it runs on another machine should 
not be the primary basis for the decision. Consideration 
must be given the machines, their I/O devices, the re- 
spective programming languages, the amount of change that 
will be required' in the code, the storage areas available, 
and the speed of execution. If the code is well documented, 
it may well be better to rewrite the code and avoid the 
pitfalls which inevitably accompany conversions of large 
complex codes. 

B. IDEAS TO PURSUE 

1 . Storage of Code on Disc 

The most efficient method of entering the code 
into the CPU is by creating a load module library which 
is discussed in reference 8. A load module is the output 
of the linkage editor. For execution, a call is made to 
this library and the code is immediately placed in core and 
is ready for execution. The only additional cards to be 
manually entered will be the input data cards. 

2 . Time Sharing 

A full discussion of time sharing is beyond the 
scope of this paper. However, its possible use should be 
mentioned at this point. Time sharing is available during 
specific periods of the day and a user may sign up for 30 
minutes at a time. During time sharing the operator may 
communicate directly with the machine. This would allow 
reinsertion of the manual control features of the code which 
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had to be removed for operation on OS. 

Space requirements would be the limiting consider- 
ation. In its present form the code is too large to fit in 
the files allocated to a particular terminal. It is possible 
that a special authorization and allocation could be ob- 
tained to allow operation on the time sharing system. 

Time sharing operation would allow the operator to 
monitor critical variables during the execution of the code 
and to observe any instabilities as they develop. This 
would allow early diagnosis of the problems associated with 
the instability and save much computing time that would be 
wasted on the OS system. 

3. Reduction of the Source Code 

Lindemuth's source code was written towards the 
goal of successful execution and had not been fully optimized 
when received at the Naval Postgraduate School. If time and 
space become critical, many of the procedures can be re- 
written in shorter form or in a time-saving form. This 
would be a time-consuming task initially, but would provide 
great savings in core storage and execution time in sub- 
sequent runs. 

Another space-saving recommendation is to specialize 
the program. In its present form the code is written to 
solve many general problems. After the code is fully de- 
bugged and operating properly, it would be advantageous 
to rewrite the code for a specific problem, for instance 
the laser produced plasma problem. This would allow all 
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of the coding which was not used for the specific problem 
to be removed and thus reduce the size of the storage area 
required for execution. 

4, Data Retention 

Presently the code stores the values of all the 
dependent variables at each mesh point for every mth time 
step as specified on input card number 3* This can 
quickly consume large areas of storage. If retention of 
this data is required for future data analysis, consideration 
should be given to storing the information on larger and 
less frequently used external storage devices. Magnetic 
tapes would be good for this purpose since the user could 
obtain his own personal tapes and use them in any way he 
required. Data could be placed on these tapes by either 
of two methods. 

If long runs are expected (for example runs of from 
one to two hours during the night or on week-ends) the 
code could be modified to write file RECORD on the tape 
instead of on the disc as is presently done. 

If short runs are expected (especially if time 
sharing can be implemented) it would be advantageous to 
write on the disc for several sessions and then use a 
utility routine to transfer the data to a tape. This is 
because the access time to the disc is much shorter than 
to the tape and this time could be used for execution during 
a short session. 
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APPENDIX A 



MACHINE DIFFERENCES 



This program was originally written in CDC FORTRAN for 
use on the CDC-3400 and later midified to run on the 
CDC-6400. While FORTRAN is generally regarded as a uni- 
versal programming language, there are a great many 
differences in the versions written for machines of differ- 
ent manufacture. The following list compares major 
differences between the CDC 3^00/6^00 and the IBM 3^0 
FORTRAN versions which were discovered during the con- 
version work on Lindemuth's code. Corrective action taken 
for each item follows the description. 

1. Identifiers 

a) CDC allows up to eight alphanumeric characters. 

b) IBM allows up to six alphanumeric characters. This 
required locating all identifiers over six characters in 
length, renaming them with six character strings, and 
entering these changes in the source deck. Table I 

lists all of the oversize identifiers and their new names. 

2. COMMON Blocks 

a) CDC allows alphanumeric identifier names up to 
eight characters in length, or numeric identifiers up to 
eight characters. 

b) IBM allows only alphanumeric identifiers up to 
six characters in length, the first character being 



24 



TABLE I 



IDENTIFIER NAME CHANGES 



ORIGINAL 




MODIFIED 


ORIGINAL 




MODIFIED 


BCRASHK 




BCRSHK 


DBPEPDZ 




BPEPDZ 


BCZASHK 




BCZSHK 


DBREPDR 




BREPDR 


BEXINIT 




BEXINT 


DBREPDZ 




BREPDZ 


BFROMSI 




BFRMSI 


DBZEPDR 




BZEPDR 


BOUNDRY 




BOUNDY 


D3ZEPDZ 




BZEPDZ 


BZEINIT 




BZEINT 


DSI EPDR 




S I EPDR 


CBERMZM 




3ERMZM 


DSIEPDZ 




SIEPDZ 


CBERPZM 


• • 


BERPZM 


DTCDVCX 




DTDVDX 


CBERMZP 


• • 


BERMZP 


GAMIOU2 




GM1U02 


CBERPZP 


• • 


BERPZP 


IBEDATA 




IBEDAT 


CLNLAM2 




CLNLM2 


INITCON 




INTCON 


CTCRMZM 




TCRMZM 


IPLAVAC 




IPLAVC 


CTCRMZP 




TCRMZP 


I ST E PNO 




I S TP NO 


CTDRMZM 




TDRMZM 


ITCDATA 




ITCDAT 


CTDRMZP 




TDRMZP 


I TC VERS 




ITCVER 


CTDRPZM 




TDRPZM 


I OUT 1MW 




IOTIMW 


CTDRPZP 




TDRPZP 


JSTEPNO 




JSTPNO 


D2BPDR2 




D2BPDR 


LNLAM2E 


• • 


LNLM2E 


D2BRDR2 




D2BRDR 


LNLAM2I 


• • 


LNLM2I 


D2BZDR2 




D2BZDR 


NOT C DRV 


• • 


NTCDRV 


D2R0DR2 




D2RODR 


OUTPUT1 




OUTPT1 


D2SIDR2 




D2SIDR 


OVER ISU 




OVRISU 


D2TEDR2 




D2TEOR 


PRINCHR 




PRNCHR 


D2TIDR2 




D2TIDR 


RORINIT 




RORINT 


D2VPDR2 




D2VPDR 


RUNDATA 




RUNDAT 


D2VRDR2 




D2VRDR 


SETKEQO 




STKEQO 


D2VZDR2 




D2VZDR 


STARTUP 




STRTUP 


DBPEPDR 




BPEPDR 


STARTYM 




STRTYM 
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alphabetic. Lindemuth's code used numbered COMMON blocks 
(COMMON/l/) • This required changing all COMMON block 
identifiers to an alphanumeric identifier. This was 
accomplished by prefixing each numeric identifier with 
a ’C* (C0MM0N/C1/). 

3. Multiple Commands 

a) CDC allows more than one command to be placed on 

one card by placing a '$• between the commands. For 
example: A=1 $ B=2 $ C=3 

b) IBM allows only one command per card. Due to the 
large number of multiple command cards, a short program 
was written to locate these cards, separate the commands, 
place them on separate cards, and renumber the source 
deck. 

4. Computed GO TO 

a) CDC version takes the form GO TO (a,b,c) X. 

b) IBM version takes the form GO TO (a,b,c), X. This 
required locating all of the computed GO TO commands and 
inserting a comma after the right par&nthesis . 

5. Floating Point Constants 

a) CDC allows floating point constants to range from 
lo + 3 22 to 10' 2 93. 

b) IBM allows floating point constants in the range 
10 ±75 . While the results of calculations are not outside 
the range, intermediate values of some calculations do 
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fall beyond this range resulting in UNDERFLOW or OVERFLOW 
and subsequent termination of program execution. This 
was corrected by locating such calculations during test 
runs and rewriting the functions such that overflow or 
underflow would not occur. 

6. I/O Commands 

a) CDC I/O commands take the form READ/WRITE INPUT 
TAPE a, b, list where *a' is the data set identifier, 'b' 
is the format statement number, and 'list* is the set of 
variables to be read or written, 

b) IBM I/O commands take the form READ/WRITE (a, b ) list 
where 'a', *b', and 'list' are defined above. This required 
locating all READ/WRITE INPUT TAPE commands and rewriting 
them in the IBM format. It was not practical to write a 
program to accomplish this because the CDC I/O command 

did not always appear in the same location on the card. 

7. Transferring Bulk Data 

a) CDC FORTRAN has a data transfer command of the 
form BUFFER IN/OUT (a,n)(A,B) where 'a' is the device 
number, 'n' is the direction mode, 'A* is the first word 
to be read or written, and 'B' is the last word to be 
read or written. This command causes all of the words 
between 'A' and *B* to be transferred into or out of core 
in the binary mode. The parameter 'n' will specify 
whether data is transferred beginning with 'A' and ending 
with 'B' or vice versa. 
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b) IBM FORTRAN provides no such command. To transfer 
binary data in or out of core is accomplished by a command 
of the form READ/WRITE (n) list where 'a' is the device 
number and 'list* is the list of words to be transferred. 
Each item must be either listed explicitly or implicitly 
within an implied DO loop. Conversion required locating 
each BUFFER IN/OUT command and changing it to the IBM 
READ/WRITS (a) command. Again it was not feasible to write 
a program for machine conversion. 

8. I/O Buffering Check 

a) CDC checks the status of a EUFFER IN/OUT operation 
with a command of the form IF (UNIT, a) n,m where 'a' is 

the device number and *n ' and 'm' are statement numbers. 
This command looks at unit 'a* which had previously been 
given a BUFFER command. If the unit is still transferring 
data, program control is transferred to statement number 
'n*. If the buffering is complete, program control is 
transferred to statement number 'm'. This allows the CPU 
to work on another section of the program which does not 
involve the data being transferred. In this way, the CPU 
is not tied up waiting for the slow i/o function to be 
completed. 

b) IBM does not have this command. Conversion con- 
sisted fo simply deleting these commands from the source 
deck. 
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9. Files 

a) CD C FORTRAN contains a subroutine which will create, 
destroy, open, or close a file during execution of the 
program. Input parameters for this subroutine may be 
variable identifiers. 

b) IBM FORTRAN has no such subroutine. File size 
and description must be determined by the operator before 
the program is compiled. File declarations must be entered 
on the JCL cards. Conversion consisted of deleting all 
calls to this subroutine from the source deck. In addi- 
tion, special instructions for creation of files and data 
sets were written. 

10. Word Length 

a) CDC machines use 48 bit words. 

b) IBM machines use 32 bit words. This has indirect 
effects which are included in other subjects below, 

11. Hollerith Assignment 

a) CDC FORTRAN allows assignment statements of the 
form IW0RD=8Hword where 'word' is a Hollerith string of 
eight characters. 

b) IBM FORTRAN does not allow direct Hollerith 
assignment. Such assignments must be made with a DATA 
statement. Since many of the assignments were to 
dimensioned variables, the following conversion was used. 

All Hollerith strings were assigned to identifiers in a 
DATA statement. The identifiers were of the form C0MT1, 
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C0MT2, and so forth. In the program the subscripted 
variables were assigned the new identifiers instead of the 
Hollerith strings. An example is: IDA( 1 )=C0MT1 where 

COMTl='word' instead of IDA( 1 )=8Hword . 

12. Hollerith Word Size 

a) Because of the 48 bit word and BCD format, CDC 
can assign an eight character Hollerith string to a 
REAL*4 integer variable. 

b) Because of the 32 bit word and EBCDIC format, IBM 
can only assign a four character Hollerith string to a 
REAL*4 variable. A two step conversion was required. 
First, all variables containing Hollerith strings were 
declared REAL* 8 . This satisfied assignment requirements. 

A second problem arose in COMMON/C 3/ • The IDA array 
is equivalenced to DA(71) in this common block. The DA 
array is type REAL*4 and the IDA array had been changed 
to type REAL*8. Therefore the DA array had to be extended 
by forty words to accommodate the double precision IDA 
array. In addition all variables which were equivalenced 
to DA(lll) or higher had to be equivalenced forty words 
higher to avoid interference with the expanded IDA array. 

13. Sense Switches 

a) CDC 3400 provides a capability of entering a 
'O' or • 1* into specified locations during execution of 
the program without a PAUSE or program stop. This is 
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accomplished by manually positioning switches at the 
operators console. The command IF SENSE SWITCH (n) m 
will cause program control to be transferred to statement 
number *m* if snese switch 'n' is on. Otherwise control 
will come from the next executable command. 

b) No sense switch capability as described above 
is available to the general user on the IBM operating 
system. Conversion consisted of removing all sense switch 
commands from the source program, along with much of the 
diagnostic and manual control capability. Presently the 
operator must know what diagnostics or special control 
features he wants and enter these requests on input data 
cards before he runs the program. 

14. Pause 

a) CDC allows use of the PAUSE command which will 
cause the program to halt execution and await instruction 
from the operator's console. This command was used ex- 
tensively with the diagnostic subroutines and manual 
control options mentioned in paragraph 13 above. 

b) IBM 360 has the PAUSE capability, but does not 
permit its use by users on the operating system (OS). 

If the code can be made compatable with the CP/CMS system 
(time sharing), the PAUSE command can be used to advantage 
as intended. 
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15* Character Representation 

a) Lindemuth’s source code was punched in BCD. 

b) Although the IBM machine can convert from BCD 

to EBCDIC for computation, all of the periferal equipment 
uses EBCDIC representation. Conversion of the source 
deck was made by the same program discussed in paragraph 3 
above . 

1 6. Comment Card Lists 

a) CDC allows an unlimited number of consecutive 
comment cards to be placed in a source program deck. 

b) IBM allows only 31 consecutive comment cards to 
be placed in a source program deck. If more than 31 
comment cards are in a program deck during compilation, 

31 cards will be listed, followed by an error message 
terminating the listing of the comment cards. Compila- 
tion will proceed with the next executable statement. 

Lists longer than 31 cards may be printed by inserting 
an executable statement every 31 cards. The command 
CONTINUE is sufficient for this purpose. 

17. Format Statements 

a) CDC FORTRAN limits a Hollerith string within a 
format statement by inclosing the string with asterisks. 

b) IBM limits the string with single quotes. Con- 
version required changing all of the limiting asterisks 
to single quotes. 
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18. Integer Size 

a) CDC’s 48 bit word allows integers up to 2^-1. 

b) IBM’s 32 bit word allows integers up to 

This caused integer constants used to set boundary condi- 
tions to be truncated and thus set improper boundary 
conditions. Conversion required changing the input format 
so that the boundary condition constant was read as two 
five-digit words instead of the previous ten-digit word. 
The scheme for extracting the boundary condition array 
values from the input constant was changed to operate on 
the two smaller words. 

19. Subroutine Names 

a) CDC allows alphanumeric identifiers up to eight 
characters in length. 

b) IBM allows alphanumeric identifiers up to six 
characters in length. This required locating all of the 
subroutines with names greater than six characters and 
rewriting the identifiers with six character strings. 

The following list contains the subroutine names which 
were too long and the new six character names which were 
assigned to them. 



BCRASHK 


changed to 


BCRSHK 


BCZASHK 


BCZSHK 


BEXINIT 




BEXINT 


BFROMSI 




BFRMSI 


BOUNDRY 




BOUNDY 


OUTPUT 1 




0UTPT1 


OVERISU 




OVRISU 


RUNDATA 




RUNDAT 


STARTUP 




STRTUP 
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APPENDIX B 



JOB CONTROL LANGUAGE (JCL) CARDS 

A. USE IN THE OPERATING SYSTEM (OS) 

Job control language (JCL) is the programming language 
which controls the IBM system J60 machine. It is the 
primary language, and other programming languages such as 
FORTRAN, ALGOL, COBOL, and so forth are used as sub- 
programs. The JCL program is concerned with setting up 
the machine to handle a specific program which is written 
in another language. This includes keeping track of account- 
ing data, calling special programs such as the FORTRAN 
compiler, defining and constructing data sets, and term- 
inating the job. A more detailed discussion of JCL is 
presented in reference 3» 

B. REQUIRED CARDS 

There are four basic types of JCL cards. The following 
are brief descriptions of their functions. A detailed 
explanation of their functions is contained in reference 3. 

1 . Job Card 

The job card is the first card of the program and 
contains all the accounting data for the particular job. 

2. EXEC Card 

This is the second card of the program and contains 
instructions for execution of the program. This card can 
be used to call a catalogued procedure or a program stored 
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in a private library for execution. It is also used to 
specify various parameter options which are available in 
catalogued procedures. REGION, or core storage, require- 
ments can also be specified on this card. 

3. Data Definition (DP) Card 

This card follows the EXEC card and is used to set 
up any data sets that may be required by the program. 

It specifies such parameters as the data set identifiers, 

device on which the data set will reside, amount of memory 
which will be required for the data set, disposition of 
the data set, and the data control block (DCB ) which 
describes the construction of the data set. 

4. Delimiter (/*) Card 

This card limits the extent of a particular operation. 

It is placed at the end of a set of cards which are being 

read for a particular operation in the program and signals 
the end of that particular set. For instance, when a 
source program is compiled, the /* card is placed at the 
end of the source deck and signals the end of the cards 
which are to be compiled. If data cards are to be input 
for execution of the source program, they are placed be- 
tween a DD card and another /* card which this time signals 
the end of the input data set. 



35 



C. CATALOGUED PROCEDURES 



1 . Descriptions and Use 

Catalogued procedures are programs which are 
written and placed in the machine library for use by the 
general user. These programs are catalogued because they 
are very common programs and are used frequently by many 
users. Each catalogued procedure is named with an identi- 
fier. To call a specific procedure, its identifier is 
placed on the EXEC card following the EXEC. For example, 
the EXEC card 
// EXEC FORTC 

would call the catalogued procedure FORTC which would then 
•compile only* a FORTRAN program source deck which followed 
in the input stream. 

2 . Over-ride Procedures 

Since the catalogued procedure is a program which 
is already written, all available options are preset. These 
may be modified by over-riding the particular preset option 
as described below. 

a. Time Over-Ride 

The general user at NPS is currently allowed 20 
seconds for his job. This includes compilation and execu- 
tion. It does not include the time used in reading or 
writing on data sets. This 20 second time may be over- 
ridden by specifying the required time on the JOB card 
following the •name* field. For example, the following 
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JOB card 

//YAG30833 JOB (0833, 0739FP.FA11), ’THESIS’ ,TIME=10 
would allocate ten minutes to job YAG30833. 

b. REGION Over-Ride 

The general used is allowed 100K bytes of 
magnetic core for his job. This includes space for the 
program and any external subroutines used by the main 
program. The REGION over-ride is made on the EXEC card. 
For example the EXEC card 

//EXEC FORTCLG , REGION. F0RT=150K 
would allocate 150K bytes of core to the FORT step in the 
catalogued procedure FORTCLG. 

c. Data Set Space Over-Ride 

The space allocation varies with the data set. 
For example the general user is allowed three cylinders 
of 2314 disc storage for data set FT06F001, which is the 
data set for the line printer. If more space is required, 
the over-ride is placed in the SPACE parameter of the DD 
statement for that data set. For example the statement 
//GO.FT06F001 DD SPAGE-(CYL, (3,1) ) 
would initially allocate three cylinders of disc storage 
to the data set. If this were not sufficient, additional 
space would be allocated one cylinder at a time until a 
total of 18 cylinders is used. 

d. Parameter (FARM) Over-Ride 

Standard options are set in a catalogued pro- 
cedure when it is written. These are outlined in a 
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description of the procedure. If any of these options 
need to be modified, they may be over-ridden on the EXEC 
card. For example the catalogued procedure FORTCLGP has 
the standard options WAP and NODECK in addition to a list 
of others. The MAP option prints out a map showing the 
location of all identifiers used in the program. NODECK 
indicates that no object deck is to be punched. If the 
operator wants the procedure to punch an object deck in 
addition to the normal procedure, the following card 
// EXEC FORTCLGP, FARM. FORT=' MAP, DECK’ 
would accomplish this operation. The standard compiler 
option is NOMAP. Procedure FORTCLGP has changed this to 
MAP. Note that it was necessary to specify MAP again in the 
over-ride statement. When one member of the parameter 
list is over-ridden , all of the other members revert to 
the standard options, so MAP would have become NOMAP if 
it had not been respecified. 

3. QUICKRUN 

QUICKRUN is a special inplementation of the 
catalogued procedure FORTCLG. It is available to general 
users whose jobs run 20 seconds or less, use 100K bytes 
of core or less, do not use special periferal equipment, 
and do not call specialized library subroutines. 

The procedure reads, compiles, and executes a FORTRAN 
source program deck. The output suppresses all page skips 
unless the page skip suppress is over-ridden. The advantage 
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of a QUICKRUN job is its fast turn-around time of about 
five minutes since these jobs have priority over all other 
classes of jobs. This mades QUICKRUN an ideal procedure 
to use for testing or debugging small programs or sub- 
routines. The JCL deck setup for a QUICKRUN job is listed 
in paragraph D of this appendix. 

D. USEFUL PROGRAMS 

The following JCL programs were used in the conversion 
of Lindemuth's code. Each program is described briefly 
and is then listed as it would appear for execution. The 
examples all use the author's accounting data. The reader 
must insert his own accounting data where appropriate. 
Reference 3 fully describes the required accounting data. 

In these programs the option 'MAP' refers to the listing 
of the contents of storage areas by identifier names. This 
option is standard on all of the catalogued procedures. 

The option 'LIST' causes the object deck listing to be 
printed out. This listing is useful in debugging because 
it associates relative storage addresses with statement 
numbers in the source program. 
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THIS PROGRAM WILL COMPILE THE SOURCE PROGRAM DECK FOLLOWING 
THE //FORT . SYS I N DO * CARD AND WILL PRODUCE AN 
OBJECT DECK. 

/ / YAG3 0833 JOB ( 08 33 , 073 9FP , F A1 1 ) , • C HARL I E YAGER THESIS* 

// EXEC FORTCD 
//FORT .SYSIN DO * 

SOURCE PROGRAM DECK. 

/* 



THIS PROGRAM WILL COMPILE THE SOURCE PROGRAM DECK FOLLOWING 
THE //FORT. SYSIN DD * CARD AND WILL ADD A PRINTED LISTING 
OF THE OBJECT DECK MAP. 

//YAG30833 JOB ( 0833 , 0739FP , FA 1 1 ) , 'CHARL I E YAGER THESIS* 

// EXEC FORTC,PARM.FCRT='MAP,LIST« 

//FORT. SYSIN DD * 

SOURCE PROGRAM DECK. 

/* 



THIS PROGRAM WILL COMPILE THE SOURCE PROGRAM DECK FOLLOWING 
THE //FORT. SYSIN DD * CARD' AND WILL PROVIDE BOTH THE 
OBJECT DECK AND ITS MAP LISTING. 

// YAG3 0833 JOB (0833,0739FP,FA11), 'CHARLIE YAGER THESIS' 

// EXEC FORTCD,PARM.FORT='MAP, LIST, DECK' 

//FORT .SYSIN DD * 

SOURCE PROGRAM DECK. 

/* 



THIS PROGRAM WILL COMPILE, LINK EDIT, AND EXECUTE THE SOURCE 
PROGRAM FOLLOWING THE //FORT. SYSIN DD * CARD. 

// YAG3 0833 JOB ( 0833 , 0739F P , F All ) , ' CHARL I E YAGER THESIS' 

// EXEC FORTCLG 
//FORT. SYSIN DD * 

SOURCE PROGRAM DECK. 



/* 

//GO. SYSIN DD * 

DATA DECK IF REQUIRED. 



/* 



THIS PROGRAM WILL COMPILE, LINK EDIT, AND EXECUTE THE SO^AC£ 
PROGRAM FOLLOWING THE //FORT. SYSIN DD * CARD AND WILL 
PUNCH AN OBJECT DECK. 

//YAG30833 JOB ( 0833 , 0739FP , FA11 ),' CHARL I E YAGER THESIS' 

// EXEC FORTCLGP,PARM.FORT='MAP,DECK' 

//FORT .SYSIN DD * 

SOURCE PROGRAM DECK. 

/* 

//GO. SYSIN DD * 

DATA DECK IF REQUIRED. 

/* 
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THIS PROGRAM WILL COMPILE, LINK EDIT, AND EXECUTE THE SO U&CE 
PROGRAM FOLLOWING THE //FORT.SYSIN DD * CARD AND WILL 
PROVIDE AN OBJECT DECK MAP LISTING. 

// YAG3 0833 JOB ( 0833 , 0739FP , FA1 1 ) , • C HARL I E YAGER THESIS* 

// EXEC FORTCLG,PARM.FORT=*MAP,LIST* 

//FORT.SYSIN DD * 

SOURCE PROGRAM DECK. 

/* 

/ /GO .SYS IN DD * 

DATA DECK IF REQUIRED. 

/* 



THIS PROGRAM WILL COMPILE, LINK EDIT, AND EXECUTE THE SOURCE 
PROGRAM FOLLOWING THE //FORT.SYSIN DD * CARD AND WILL 
PROVIDE BOTH THE OBJECT DECK AND ITS MAP LISTING. 

//YAG30833 JOB ( 0833 , 0739FP , FA11 ) , 'CHARLIE YAGER THESIS' 

// EXEC FORTCLGP, P ARM. FORT= ' MAP , LIST , DECK* 

//FORT.SYSIN DD * 

SOURCE PROGRAM DECK. 

/* 

//GO .SYS IN DD * 

DATA DECK IF REQUIRED. 

/* 



THIS PROGRAM WILL LINK EDIT, AND EXECUTE THE SET OF OBJECT 
DECKS FOLLOWING THE //LINK.SYSIN DD * CARD. 

//YAG30833 JOB ( 0833 , 0739F P , FA1 1 ) , ' CHARL I E YAGER THESIS' 

// EXEC FORT LG 
//LINK.SYSIN DD * 

SET OF OBJECT DECKS. 

/* 

//GO. SYS I N DD * 

DATA DECK IF REQUIRED. 

/* 



THE FOLLOWING PROGRAM WILL COMPILE, LINK EDIT AND EXECUTE 
THE SOURCE PROGRAM FOLLOWING THE JOB CARD USING THE QUICKRu^ 
FACILITY. NOTE THAT THE JOB CARD IS OF A DIFFERENT FORMAT 

// QUICKRUN 

#YAG3083 3 FORTGO ( 0833 , 0739FP , FA1 1 ) , • QUI CKRUN EXAMPLE* 

SOURCE PROGRAM DECK. 

# DATA 

DATA DECK IF REQUIRED. 

// 
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THIS PROGRAM WILL LIST THE CARDS FOLLOWING 
THE //SYSUT1 DD DAT A, DCB=BLKS IZE = 80 CARD. 

//YAG30833 JOB ( 0833 , 0739FP , FA 11 ) t 'CHARL I E YAGER THESIS' 
// EXEC PGM= I E BG EN ER 
//SYSPRINT DD SYSOUT=A 
//SYSIN DD DUMMY 

//SYSUT2 DD SYS0UT=A,DCB=(RECFM=F,BLKSIZE=80) 

//SYSUT1 DD DATA»DCB=BLKSIZE=80 
INPUT DATA DECK. 

/* 



THIS PROGRAM WILL COMPILE THE SOURCE PROGRAM DECK FOLLOWING 
THE //FORT. SYSIN DD * CARD. 

//YAG30833 JOB ( 0833 , 0739FP , FA11 ),' CHARL I E YAGER THESIS' 

// EXEC FORTC 
//FORT. SYSIN DD * 

SOURCE PROGRAM DECK. 

/* 
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APPENDIX C 



FILES 



A. DEFINITIONS 

A distinction should be made between FILES and DATA 
SETS. This paper will use the following descriptions? 

1. File 

A file is a collection of data or information which 
is used for a specific purpose. This collection will be 
given a characteristic identifier name which can be used 
for reference. For example, the file to hold all the data 
generated by the program is given the identifier DATA. 

2 . Data Set 

A data set is the device from which or to which 
data and information is passed. Data sets are defined in 
JCL cards which are not part of the source program code. 
Common data sets are the card reader, card punch, line 
printer, magnetic tapes, magnetic discs, and paper tapes. 
See appendix H for a detailed description of data sets 
and their use. 

B. DESCRIPTIONS 

1. DATA 

This file consists of all the input data cards used 
to start or restart a program. In its present form the 
program will use this file with data set FT05F001, the card 
reader. Appendix E describes all of the input data cards 
used by this program. 
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2. MHDOUT 



This file consists of all the output information 
generated during execution of the program. In its present 
form this program will use this file with data set FT06F001, 
the line printer. 

3. RECORD 

This file consists of all the information necessary 
to restart a program plus the values of each variable at 
each mesh point at each time step. Figure 1 is a graphical 
description of the organization of the file. In its pres- 
ent form this program will use this file with data set 
FT10F001 , a data set on disc storage. 



Begin 



W 

R 

I 

T 

E 



End 

Figure 1. Organization of File DATA 



End of File 


175 WORDS 


Record Gao 


NDR+NDZ 

WORDS 


Record Gao 


25+NDR*NDZ*NV 

WORDS 


Record Gao 


25+NDR*NDZ*NV 

WORDS 


Record Gap 


• 

• 

• 


Record Gao 


25+NDR*NDZ*NV 

WORDS 


End of File 



All 175 words of DA array 
Only words used of RZ array 



Initial DADT and VAR arrays 
All 25 words of DADT 
Only words used in VAR 



DADT and VAR after 1 time step 
Many time steps 

DADT and VAR after ' NDT ' time steps 
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k. OLDREC 



This file is identical to RECORD. OLDREC is used 
in the restart mode of operation and is created from RECORD 
by changing the name of RECORD to OLDREC. When the program 
is operated in the restart mode, OLDREC is used by sub- 
routine STRTUP as a source of input (note that this is 
not a part of file DATA) and the program will enter the 
new computed values into a new file RECORD. In its pres- 
ent form this program will use this file with data set 
FT30F001, a data set on disc storage. 

5. SCRAT 

This file is a work area for subroutines CHECK 
and SPECHK and is created only when these diagnostic 
routines are used. In its present form this program will 
use this file with data set FT40F001, a data set on disc 
storage which is destroyed after each run of the program. 

C. USE WITH READ/WRITE COMMANDS 

Data sets are written onto and read from their re- 
spective units by the standard FORTRAN statements READ and 
WRITE. Information buffered in and out may take two forms, 
either format controlled or binary. 

The format controlled READ statement takes the form 
READ(a,b) list, where 'a* is the device number (which is 
identical to the integer in the third and forth positions 
of the data set number), 'b* is the format statement number, 
and 'list* is the set of variables to be read. 
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The binary READ statement which has no format control, 
takes the form READ (a ) list where 'a* is the device number 
and ’list' is the set of variables to be read. The binary 
READ/WRITE is used to transfer information between the data 
set and core without changing the bit pattern. 

Files DATA and MHDOUT use format controlled I/O state- 
ments while files RECORD, OLDREC , and SCRAT use binary 
I/O statements. 

This program uses integer identifiers for the device 
number in all of the I/O statements. The identifiers are 
set to the appropriate integer values in the initial steps 
of the main program. This allows flexibility in the choice 
of data sets in that a data set can be altered by changing 
one card in the main program instead of changing every 
READ/WRITE card in the entire program. Table II is a ref- 
erence of files and data sets including all identifiers 
which are used by the program in its present form. 



46 



DATA SET IDENTIFICATION REFERENCE 
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APPENDIX D 



SUBROUTINE DESCRIPTIONS 

This section will list all of the subroutines included 
in the present code and will briefly describe their 
functions and operations. The subroutine functions can be 
divided into five categories. 

1. LOGIC AND CONTROL : These are the subroutines which 

primarily control the operation of the code by directing 
the call of other subroutines after manipulation of control 
variables. Subroutines included in this category are MAIN, 
BCRETS, BCZOTS , BCRSHK , BCZSHK, BOUNDY, and OVRISU. 

2. INITIALIZATION : These are the subroutines which 

set up all of the initial values required to run a particular 
problem. Except for three logic control cards, all of the 
input data cards are read by the subroutines in this 
category. Subroutines included in this category are BEXT , 
INIT1, INIT4, INIT5, MATRIX, MESH, SETUP, STRTUP, and 
TCINIT. 

3. OPERATION: These are the subroutines which perform 

the actual calculations during execution. Subroutines in- 
cluded in this category are BCR, BCZ, BFRMSI , M.AT2 , ONESID, 
SHOCK, TRANCO, and TRIANG. 

4. OUTPUT : These are the subroutines which process 

information for inclusion in the output files RECORD and 
MHDOUT. While much of their processing is done before 
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any calculation takes place, their operations are still 
considered output. Subroutines included in this category 
are BEXINT, BUFFER, 0UTPT1 , PLOTR , RUNDAT , TEPLOT, and 
ZBC. 

5. DIAGNOSTIC s These subroutines are used to diagnose 
numerical operations. By performing inverse operations 
on calculated matrices, variations in the bit patterns of 
words can be observed and analyzed. The subroutines in 
this category are CHECK and SPECHK. 

The following list provides a brief description of the 
function of each subroutine. 

A. MAIN PROGRAM 

This is one of the logic and control routines and 
controls the general operation of the entire program. It 
reads cards number 1, 2, and 3 from the input data set 
to control diagnostics requests and data retention in file 
RECORD. This routine also checks for negative density, 
ion temperature and electron temperature. If these 
quantities occur, the program will terminate and print out 
that negative values were found during execution in a par- 
ticular time step, 

B. BCR (Boundary Conditions for R) 

This is one of the operation subroutines and is 
called by MAIN and BOUNDY. BCR sets up boundary matrices 
and vectors along the and R NDR lines corresponding to 
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equations 3*5-3 and 3*5-^ of reference 1. Boundary 
condition control parameters are passed from the initializing 
subroutines through COMMON blocks C4 and C 6 and determine 
which of a set of pre-defined boundary conditions will be 
applied to a particular variable. Subroutine RBC forms 
printed statements of the selected boundary conditions. 

C. BCRETS (Boundary Conditions for R, Even Time Step) 

This is one of the control and logic subroutines and 

is called by MAIN. It calls RBC through BOUNDY and is con- 
cerned with boundary calculations corresponding to equations 
3*5-3 an d 3*5-^ of reference 1 with 'n* replaced by *n+l' 
and 'n+1' replaced by *n+2*. Again operations are on lines 

R 1 and R NDR’ 

D. BCRSHK (Boundary Conditions for R, after SHocK 

SmoothingJ 

This is one of the control and logic subroutines and 
is called by MAIN, It calls RBC through BOUNDY and is con- 
cerned with boundary calculations along line R^ and R^^ * 
BCRSHK causes boundary values to be readjusted after shock 
smoothing if used. 

E. BCZ (Boundary Conditions for Z) 

This is one of the operation subroutines and is 
called by MAIN and BOUNDY. BCZ sets up boundary matrices 
and vectors along the Z l and Z NDZ lines corresponding to 
equations 3*5-5 and 3*5-6 of reference 1. Boundary 
condition control parameters are passed from the initializing 
subroutines through COMMON blocks C4 and C6 and determine 
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which of a set of pre-defined boundary conditions will be 
applied to a particular variable. Subroutine ZBC forms 
printed statements of the selected boundary conditions. 

F. BCZOTS (Boundary Conditions for Z, Odd Time Step) 
This is one of the control and logic subroutines and 

is called by MAIN. It calls ZBC through BOUNDY and is con- 
cerned with boundary calculations corresponding to equations 
3.5-5 and 3.5-6 of reference 1 with 'n* replaced by 'n+1' 
and 'n+1' replaced by 'n+2'. Again operations are on lines 

Z 1 and Z NDZ’ 

G. BCZSHK (Boundary Conditions for Z, after SHocK 

Smoothing! 

This is one of the control and logic subroutines and 
is called by MAIN. It calls ZBC through BOUNDY and is con- 
cerned with boundary calculations along lines Z 1 and % m z' 
BCZSHK causes boundary values to be readjusted after shock 
smoothing if used. 

H. BSXINT(J»K) (B Field External, INiTially) 

This is one of the output subroutines and forms 

the labels for the external B field values which are read 
in on input card number 10. Data values are passed to other 
subroutines through COMMON/C 14/. 

I. BEX T (B Field EXTernal) 

This is one of the initializing subroutines and 
works in conjunction with BEXINT to provide a uniform ex- 
ternal B field. BEXT- is used when IBE/0. See section 3.8, 
page 83 of reference 1. Initial values are passed to BEXT 
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through COMMON/C 14/ and the calculated values are returned 
to other subroutines through COMMON blocks Oil and C26. 

J. BFRMSI ( J ,K , BR , BZ ) (B Field FRoM PSI) 

This is one of the operation subroutines and calcu- 
lates the B field components BR and BZ from the stream 
function . BFRMSI is called from SETUP and is used only 
when anisotropic thermal conductivities are used (see 
page 29 of reference 1 ). 

K. BOUNDY(l£NT t NNNN NNN, ID) ( BOUND aRY) 

This is one of the control and logic subroutines 
and works in conjunction with BCRETS, BCRSHK, BCZOTS, and 
BCZSHK to perform boundary calculations corresponding to 
equations 3 . 5-5 and 3 * 5-6 of reference 1 . 

L. BUFFER ( IROR , I ARRAY , ITP ) 

This is one of the output subroutines and controls 
the transfer of data into and out of file RECORD. The 
data set for file RECORD is defined for logical records 
of 100 bytes which corresponds to a list of 25 words. The 
arrays DA, RZ, and VAR are all longer than 25 words and 
require additional control to be read into file RECORD. 

This is the function of BUFFER. 

The argument IROR for READ or RITE) specifies 
whether the subroutine is to READ from or WRITE into the 
data set. If IR0R=0, BUFFER READS. If IR0R=1, BUFFER 
WRITES. 
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IARRAY specifies which of three arrays are to be 
used with the I/O command. IARRAY=1 specifies array DA, 
IARRAY=2 specifies array RZ, and IARRAY=3 specifies array 
VAR. 

ITP is the device identifier for the data set and 
is included in the argument list as the variable ITP. Its 
value is specified elsewhere in the program. 

M. CHECK 

This is one of the diagnostic subroutines and works 
in conjunction with SFECHK to print internal diagnostics. 
These diagnostics are used to determine which terms in 
the difference equations are responsible for non-physical 
values. Use of CHECK and SPECHK is determined by input 
parameters on cards 1, 2, and 13 of the input data set. 
Examples of the output from CHECK and SPECHK are included 
in reference 9. 

N. INIT1 ( iNIT ializing Routine Number 1_) 

This is one of the initializing subroutines and 
sets up the VAR array at time t=0. This subroutine sets up 
a uniform plasma such as that used in theta-pinch studies. 
This problem is discussed in chapter 5 of reference 1. 

O. INIT4 ( INIT ializing Routine Number 4) 

This is one of the initializing subroutines and 
sets up the VAR array at time t=0. This subroutine sets up 
a high density center plasma surrounded by a low density 
background such as that used in laser produced plasma 
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studies. This problem is discussed in chapter 6 of 
reference 1. Initial condition parameters are read on 
card number 12A of the input data set. 

P. INIT5 ( INIT ializing Routine Number £) 

This is one of the initializing subroutines and 
sets up the VAR array at time t=0. This subroutine sets up 
a discontinuous plasma such as that used in hydrocode 
test problems. This problem is discussed in chapter 4 
of reference 1. 

Q. MAT2 ( J ,K) 

This is the primary operation subroutine of the 
program. It determines matrices A, B, and C and the vector 
V which appear in difference equations 3.2-1, 3«5-l» and 
3.5-2 of reference 1. The matrices A’, B', and C* and 
vector V* (3*1-13» 3*3-12, and 3*3-25 of reference 1) are 
calculated and then converted to A, B, G, and V. The code 
is set up as if all ten components of W are to be calculated. 
However, if IV(K) is greater than NV, then the Kth component 
of W is not being calculated. MAT 2 is called by MIN and 
passes information back and forth through numerous COMMON 
blocks as described in appendix F. 

R. MATRIX 

This is one of the initializing subroutines and is 
the first one called by subroutine STRTUP in the start-up 
mode. It reads cards number 6, 7» and 8 from the input 
data set. The important control parameters in this 
subroutine are JAV and JVERS. 



1. JAV : JAV determines which alternate variables 

are to be used, if any. An alternate variable is one which 
can be determined from a component of the vector W by either 
multiplying or dividing by the density. JAV is the param- 
eter in a computed GO TO statement which branches to 
different values for the integer constant IIAV, IIAV 

is later used to set the values in the vector IAV which is 
then passed to other subroutines through C0MM0N/C6/. 

Figure 1 lists the values available for JAV and the function 
performed for each. 

2. JVERS : JVERS determines a preselected set of 

program options. Included in each set is the number of 
dimensions, the number of variables, the variables to be 
calculated and their order in the vector U, and the coord- 
inate system to be used in the calculation. Figure 1 lists 
the values available for JVERS and the function performed 
for each. 

3. Other Parameters: The integer constants JJDDR 

and JJDDZ are read and used to assign values to the JDDR 
and JDDZ arrays which are then passed to other subroutines 
through C0MM0N/C22/. Appendix E lists the input parameters 
read by MATRIX. 

MATRIX also forms messages in the IDA and ADDA vectors 
which are printed by RUNDAT. These messages list the input 
parameters read by MATRIX. 
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JAV 


IIAV 


Alternate Variables 


1 


0 


No Alternate Variables 


2 


1000110000 


SV=1/R0 . . . PI=TI*R0,PE=TE*R0 .... 


3 


1 




4 


110000 


. . . . PI=TI*R0,PE=TE*R0 .... 


5 


2 




6 


110002 


.... PI=TI*R0,PE=TE*R0 . . . AP 


7 . 


110001 


.... PI=TI*R0,PE=TE»R0 . . . SB=SI/R0 



JVERS 


DIM 


IC00RD 


U Vector 


0 


2 


0 


(RO ,VR ,TI ,TE ,3Z) 


1 


2 


0 


(RO,VR,TI,TE,BP,BZ) 


2 


1 


0 


(R0,VR ,TI ,TE ,BZ) 


3 


1 


0 


(RO.VR.TI ,TE,BP,EZ) 


4 


1 


0 


(RO , VR ,TI ,TE ) 


5 


2 


0 


(RO ,VR , VP ,VZ ,TI ,TE ,BP ,SI ) 


6 


2 


0 


(RO,VR,VZ,TI,TE,SI) 


7 


1 


0 


(RO , VR, TI ,TE ,BP ,SI ) 


8 


1 


0 


(RO.VR.TI, TE.SI) 


9 


2 


0 


( RO , VR , VZ ,TI ,TE ) 


10 


1 


1 


(RO.VR.TI ,TE ) 


li 


1 


-1 


(RO.VR.TI ,TE ) 


12 


1 


-1 


(RO.VR.TE) 


13 


1 


1 


(RO ,VR ,TE ) 


14 


1 


0 


(RO.VR.TE) 



Figure 1. MATRIX Control Parameter List 
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Figure 1 is to be used as a reference for the 
options which are presently available in subroutine 
MATRIX. The first block indicates what value of IIAV 
is assigned by the code for each value of JAV read in 
on input card number 7. For instance if JAV=2 is read 
in on card number 7, the code will set IIAV=100011000. 

This indicates that the first, fifth, and sixth components 
of the vector W will be treated as alternate variables. 

The effect of this is shown in the last column. The 
alternate variables are indicated and the ’ . ’ indicates 
that the variable is unchanged. Thus if JAV=2, the vector 
W would become W= (SV , VR ,VP ,VZ ,PI ,PE ,BR,BP ,BZ ,SI ) instead 
of W=s (RO, VR, VP, VZ, TI.TE.BR, BP, BZ, SI) . 

The second block indicates which dependent 
variables w ill be calculated for each value of JVERS 
read in on input card number 7. For instance if JVERS=6, 
the code would calculate (in two dimensions and cylindrical 
coordinates) the dependent variables RO ,VR ,VZ ,TI ,TE , 
and SI. The other components of W would not be calculated. 
This assumes that JAV=0. However, if JAV=2 while JVERS=6, 
the dependent variables which would be calculated are 
SV, VR, VZ, PI, PE, and SI. 
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S. MESH 



This is one of the initializing subroutines and 
sets up the mesh array RZ. Provisions are made for three 
different meshes in each dimension. Each mesh can be either 
uniformly spaced, or spaced geometrically. The code is 
limited to rectangular meshes set up in a rectangular 
domain. The parameters on card(s) 11A (and) 11B determine 
the mesh to be used. 

T. ONESID (ONE-SIDed Derivatives) 

This is one of the operation subroutines and changes 
standard centered first derivatives (equation 3.1-20 of 
reference 1) to one-sided derivatives (equations 3.8-7 and 
3.8-19). The dependent variables which are to use the 
one-sided derivatives are determined by the JDDR and JDDZ 
arrays which themselves were determined by the input 
parameters JJDDR and JJDDZ respectively. 

U. QUTPTl ( OUTP UT Routine Number 1) 

This is the primary output subroutine in the program. 
It prints out the values of each dependent variable being 
calculated at each mesh point for each time step as de- 
termined by NDTPNT. Plots of the dependent variables are 
presented as well as the axial and radial distances at 
each mesh point. 0UTPT1 is called by MAIN through OVRISU 
to print out when required. 

V. OVRISU 

This is one of the logic and control subroutines 
and serves merely to call other subroutines. It replaces 
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some of the overlays which were used in the original 
CDC 3400 version of the code. 

W. PLOTR 

This is one of the output subroutines and is used 
in conjunction with TEPLOT to plot the Z line values of 
TE (electron temperature) vs. R. Its coding is similar 
to the plotting coding used in 0UTPT1. 

X. RBC (R Boundary Conditions) 

This is one of the output subroutines and forms a 
part of the boundary condition messages which appear on the 
first page of the output printout. RBC is called by MAIN 
through OVRISU. 

Boundary condition information is passed to RBC 
through vectors IBCRR and IBCR in C0MM0N/C4/. Vector IIV in 
COMMON/C 1/ provides information concerning the variables to 
be calculated and the order in which they appear in the 
vector U. These control parameters determine the message 
which is buffered out to the output data set. 

This subroutine has no effect on the dependent 
variables and could be eliminated if hard copy of the 
boundary conditions is not required. 

Y. RUNDAT (RUN DATa) 

This is one of the output subroutines and produces 
the listing of initial conditions and boundary conditions 
to be listed on the first page of the output printout. 

The data which is printed out is stored in arrays IDA and 
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ADDA by many other subroutines during initialization 
procedures . 

Z. SETUP (J,K) 

This is one of the initialization subroutines and 
sets up the arrays in COMMON blocks C8, C9 # CIO, Cll, C15, 

C25» C26, and C27 which are used by MAT 2 . Transport co- 
efficients are calculated for the mesh point NJ.KE and its 
eight neighbors by calling TRANCO. Storage limitations re- 
quired computing the transport coefficients three times during 
each time step. On machines with sufficient core, the 
transport coefficients at all mesh points could be cal- 
culated prior to beginning a time step and thus save much 
computation time in this subroutine. 

AA. SHKMAT 

At the present time this is a dummy subroutine. 

BB. SHOCK (JJJ.KKK) 

This is one of the operation subroutines and does 
shock smoothing as discussed in section 3.6 of reference 1, 

It is called by MAIN and passes its data in and out through 

common/c 7 /. 

CC. SPECHK (SPEcial CHecK) 

This is one of the diagnostic subroutines and works 
in conjunction with CHECK to print internal diagnostics. 

These diagnostics are used to determine which terms in 
the difference equations are responsible for non-physical 
values. Use of CHECK and SPECHK is determined by input 
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parameters on cards 1, 2, and 13 of the input data set. 
Examples of the output from CHECK and SPECHK are included 
in reference 1, 

DD. STRTUP ( ISUB , IPU . JJ , KK ) (STaRTUP) 

This is one of the initializing subroutines and is 
the first one called by subroutine OVRISU in the startup 
mode. It reads cards number 4 and 5A or 5B from the input 
data set. All initialization is controlled by STRTUP as 
determined by the input cards. The argument list is used 
to pass values back to the calling program. 

STRTUP can be used to start a new run or to restart 
a previous run which had been terminated. 

1. Startup Mode: In this mode STRTUP reads cards 

number 4 and 5A. It then fills the boundary condition and 
ISHK vectors as specified by input parameters on card 5A, 
calls subroutines MATRIX, TCINIT, and BEXINT (if specified), 
and re terns the values of NDTPNT , IPU, 1, and NDTPNT-1 

to the calling program through the argument list. 

2. Restart Mode: The code provides two restart 

capabilities. The first will pick up execution at the point 
where the previous execution was terminated manually before 
NDT steps were completed. This type of termination requires 
the use of sense switch number one. Since the sense switch 
facility is not available, this restart mode will not be 
used at this time. 
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% The second restart mode picks up execution at 



any previous time step in any previous file. The file 
from which the execution is to continue must be placed in 
data set FT30F001 before the program is executed. Card 
number 5A need only contain values for IPU, I FROM , and 1ST. 
Setting I FROM not equal to zero will cause the code to 
read from data set FT30F001 the vectors DA, RZ, and the 
DADT and VAR set which correspond to the time step specified 
by 1ST. Card number 5B will now be read, and the values can 
be changed at this point since this mode always starts 
a new file. If IPU equals zero, MATRIX is called and new 
parameters may be specified via card number 6. Otherwise 
an immediate return to OVRISU will be made and the program 
will continue with the same parameters that were used in 
the previous run. 

EE. TCINIT( ITCDAT ) (Transport Coefficient INIT ialization ) 
This is one of the initializing subroutines and is 
the second one called by subroutine STRTU? in the startup 
mode. It reads card number 9 from the input data set. The 
primary control parameter is ITCDAT which is also the 
single argument in the subroutine parameter list. 

When ITCDAT is equal to zero, (the start-up mode), 
data card number 9 is read. If ITCDAT is not equal to zero, 

(a restart mode) the data card is not read. Instead the 
relevant information is extracted from the ADDA vector 
and assigned to the proper identifiers. Execution then 
continues as in the start-up mode. 
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TCINIT evaluates some fundamental constants and the 
constants which are used in the transport coefficients 
(equation 2.3 of reference 1). 

This subroutine also forms messages in the IDA and 
ADDA vectors which are printed by RUNDAT. These messages 
list the input parameters read by TCINIT and some of the 
values computed during execution. 

FF. TEPLOT 

This is one of the output subroutines and is used 
in conjunction with PLOTR to plot the Z line values of 
TE (electron temperature) vs. R. It is used to manually 
monitor TE for oscillations. If required, this can be 
easily modified to include all the dependent variables. 

GG. TRANCO(VARI ,TC ,DTCDY t IENT ,B5 t TCP) (TRANsport 

Coefficients) 

This is one of the operation subroutines and com- 
putes the values of the transport coefficients which are 
discussed in section 2.3 of reference 1. Constants in the 
transport coefficients are passed to TRANCO from TCINIT 
through COMMON/C 3/ » and the alternate variable designators 
are passed through C0MM0N/C6/. 

Dependent variable values are passed from the calling 
subroutine through the VARI vector in the subroutine ar- 
gument list. Computed transport coefficients are passed 
back to the calling routine through the vectors TC , TCD, 
and the array DTCDV in the subroutine argument list. 

External B field values are also passed to TRANCO through 
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the subroutine argument list in the BE vector. The other 
identifier in the argument list, IE NT , is a control ident- 
ifier which, if set to one, skips the coding for transport 
coefficient derivatives and leaves the DTCDV array equal 
to zero. 

HH. TRIANG(IFI) 

This is one of the operation subroutines and solves 
the equations A U=V, and A M=B by Gaussian elimination. 

The subroutine argument IFu determines which equation is 
to be solved. If IF1-1, then only the first equation is 
solved. If IF1 is not equal to one, then both equations 
are solved. 

The arrays A and B and the vector V are passed to 
TRIANG through COMMON/C?/ and the control parameter NV is 
passed through the DA vector of C0IVIM0N/C3/. 

The array M and the vector U are passed back to the 
calling subroutine through C0MM0N/C7/. Values of U are 
returned in V and values of M are returned in B. 

These arrays are discussed in section 3 of 
reference 1. 

II. ZBC (Z Boundary 2onditions) 

This is one of the output subroutines and forms a 
part of the boundary condition messages which appear on the 
first page of the output printout. ZBC is called by MAIN 
through OVRISU. 

Boundary condition information is passed to ZBC 
through vectors IBCZZ and IBCZ in C0MM0N/C4/. Vector IIV 
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in COMMON/C 1/ provides information concerning the variable 
to be calculated and the order in which it appears in the 
vector U. These control parameters determine the message 
which is buffered out to the output data set. 

This subroutine has no effect on the dependent 
variables and could be eliminated if hard copy of the 
boundary conditions is not required. 
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Figure Z> Cross-reference of Subroutines 

To find all of the subroutines called by a particular 
subroutine, locate that subroutine in the left-hand column 
and read across. 

To find all of the subroutines which call a particular 
subroutine, locate that subroutine in the top row and read 
down. 
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APPENDIX E 



INPUT DATA CARD FORMATS 



The following list includes all input cards which 
may be used in the program. They are not all required 
for any one run. The numbers assigned to these data 
cards are referenced in the subroutine descriptions of 
appendix D. 

The description of each card includes the card number, 
statement number of the appropriate READ command, and the 
card format description. A typical statement number may 
be of the form 107-2. This indicates that the READ state- 
ment is not numbered and that it is located two cards in 
front of statement number 10? in the source code. Each 
line of the format description contains the card columns 
for a particular identifier, the identifier name, the 
identifier READ format, and a brief description of the 
identifier. More detailed descriptions of the identifiers 
are given in appendix I. 

With a few exceptions discussed below, all of the input 
cards listed must be included in the input stream in the 
listed order, even though the parameters may all be zero. 
The following exceptions must be observed. 

1. If the value of NREG on card number 11A is less 
than or equal to • 1', card 11B must not be included in the 
input stream. 
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2. Cards 12A and 12B must not be read together. 

These cards are input to specific initializing subroutines 
and only one of these subroutines is used during any one 
run. Therefore the input card for the subroutine which 

is not being used must be removed from the input stream. 

3. Card number 13 may be many cards with the same 
format. This card specifies what action is to be taken 
v/hen a diagnostic call is made. There must be one of these 
'instruction' cards for each mesh point that is specified 
on card number 2. As a check, the number of input cards 
number 13 will be the same as N2BCHK which is read on 

card number 1. If N2BCHK=0, then no cards number 13 are 
required in the input stream. 
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APPENDIX F 



COMMON BLOCK STORAGE AND EQUIVALENCE 
A. COMMON STORAGE 

COMMON block storage is used to pass information between 
programs or subprograms without using a parameter list. 

This reduces the core storage required for the program 
because subprograms using the same named COMMON block use 
the same physical storage area in core. 

In IBM FORTRAN, we are allowed one un-name d COMMON 
block and as many named blocks as we need. Variables or 
arrays declared in an un-name d COMMON block may not be 
assigned values in a DATA statement. However, variables 
or arrays in named COMMON blocks may be initialized in a 
DATA statement. 

The named COMMON blocks are labeled by an alphanumeric 
identifier which may be up to six characters in length, the 
first of which must be an alphabetic character. The length 
of the COMMON block is determined by the first COMMON declar- 
ation. All COMMON declarations in subsequent subprograms 
must be less than or equal in length to the first declaration. 

Because two or more subprograms use the same storage 
area, all the variables used and stored by one subprogram 
may be used by any other subprogram which has declared 
the same COMMON block. However, care must be taken to 
insure that the variables in the COMMON declaration list 
are in the same order as they are in other lists. For 
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example, if subprograms SI and S2 both contain the state- 
ment •COMMON/Cl/ A ,B , 0 * then the same identifier in each 
program would have the same value. However, if subprogram 
S2 contained the statement • COMMON/C 1/ A,C,B* then the 
value of A would be the same in both subprograms, B in SI 
would be the same as C in S2, and C in SI would be the 
same as B is S2. 

B. EQUIVALENCE 

EQUIVALENCE statements are used to assign variables 
or arrays in the same subprogram to the same physical 
area in storage. Again this reduces the core storage 
required for the program because many variables may occupy 
the same area in core storage. Variables in COMMON may 
not be equivalenced to each other. However, any number of 
variables not in COMMON storage may be equivalenced. As 
an example, consider a program which contains the following 
statements : 

COMMON/C 3/ A(5) 

EQUIVALENCE (A(l) ,X) , (A(2) ,Y(2) ) , (A(4) ,2(1)) 

DIMENSION Y(3).Z(3) 

As a result of this EQUIVALENCE statement, the following 

set will occupy the same storage location in core: 

A( 1) and Y(l) and X 
A(2) and Y(2) 

A(3) and Y(3) 

A(4) and Z ( 1 ) 

A(5) and 2(2) 

2(3) will be assigned a storage location of its own. Any 
time an operation midifies one member of an equivalenced 
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set, the other members will automatically be set to that 
value. For example, if the program set X=l, then A(l) and 
Y(l) would also contain the value 1. 

C. USE 

Lindemuth's program makes extensive use of COMMON blocks 
and EQUIVALENCE statements. Logical and control data is 
passed between subprograms via COMMON blocks. Some COMMON 
block lists are composed of identifiers, while others 
consist of large arrays which are then equivalenced to 
individual identifiers. The arrays are used to facilitate 
rapid buffering to the I/O devices. The data transferred 
through the arrays DA, RZ, DADT, and VAR are also placed on 
permanent storage devices for later use in restarting or 
data analysis. 

D. REFERENCES 

Figure 1 is a chart which cross-references subroutines 
with COMMON blocks used in the program. By reading across 
for a particular subroutine, the COMMON blocks used in that 
subroutine may be determined. Similarly by reading down 
for a particular COMMON block, the subroutines to which 
it is common may be determined. 

Table 1H lists the COMMON blocks used in this program, 
the dimension of the block, and the variable list gener- 
ally associated with the COMMON block. A notable exception 
to this list concerns subroutine MAT 2 . In COMMON blocks 
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C8 P C9, Cll, C15, and C23 the lists consist of single 
variables instead of the vectors listed in the other 
subroutines. These variables may be cross-referenced by 
referring to the appropriate COMMON declaration in the 
source listing for subroutine MAT 2 which may be found in 
the program listing at the end of this paper. 

Table IV lists the members of the DA and DADT vectors 
and the variables to which they are equivalenced . 
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Figure 1. Cross-Reference of Subroutines 
and COMMON Blocks 

To find all of the COMMON blocks used by a particular 
subroutine, locate that subroutine in the left-hand column 
and read across. 

To find all of the subroutines which use a particular 
COMMON clock, locate that COMMON block in the top row and 
read down. 
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TABLE I II 



COMMON BLOCK CONSTRUCTION 



Cl: 



DI MENS I ON ( 25) 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



CFEQ1 

CFEQ2 

CRES1 

CRES2 

CKI1 

CK 12 

CKE2 

ITCVER 

CLNLM2 

CKIP1 



14 

15 

16 

17 

18 

19 

20 
21 
22 
23 



UO 

GAM1 

GAM2 

PI 

TCRO 

BLANK 

BLANK 

BLANK 

BLANK 

BLANK 





11 


CKEP1 


24 


BLANK 




12 


EMASS 


25 


BLANK 




13 


PMASS 






C2 : 


DIMENSION t 5) 








1 


J CHK ( 1 ) 


4 


JCHK(4) 




2 


JCHK ( 2 ) 


5 


JCHK ( 5 ) 




3 


J C HK ( 3 ) 






C3 : 


DIMENSION (2 00 +VAR(X) ) 








1-175 


DA ( 175 ) 








176-200 


DADT (25) 








201- X 


VAR ( X) 


X ( MAX ) = 


1400 


C4: 


DIMEN SI ON ( 55) 








1- 10 


IBCRR( 10) 


41- 50 


I SHK ( 10) 




11- 20 


I BCR ( 10 ) 


51- 53 


NZP ( 3 ) 




21- 30 


IBCZZC10) 


54 


IMPSHK 




31- 40 


IBCZ(IO) 


55 


JSTEP ' 


C5 : 


DIMENSION 


!( 162) 








1-162 


G ( 9 , 9 , 2 ) 






C6 : 


D I MENS I ON (31) 








1 


ISTEP 


12- 21 


IVV( 10) 




2- 11 


I V ( 10) 


22- 31 


I A V ( 10) 



C7: DIMENSION (320) 

1-100 A(10,10) 
101-200 6 ( 10 , 10 ) 
201-300 C(10,10) 



301-310 

311-320 



V( 10) 
COL ( 10) 



C8: 



DI MENSI ON ( 196 ) 



1- 4 


TCZM ( 4 ) 


29- 68 


DTCD VM ( 4 


,10) 


5- 8 


T CRM ( 4 ) 


69-108 


DTCD V ( 4 , 


10) 


9- 12 


TC ( 4 ) 


109-148 


DT CDVP ( 4 


, 10) 


13- 16 


TCRP (4) 


149-188 


DTDVDX ( 4 


,40) 


17- 20 


TCZP ( 4 ) 


189-192 


W( 4) 




21- 24 


DTCDR ( 4) 


193-196 


WW (4 ) 




25- 28 


DTCDZ ( 4 ) 
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C9 : DIMENSIONED) 





1- 10 


VARIZM( 10) 


51- 60 


DVDR( 10) 




11- 20 


VARIRM( 10) 


61- 70 


DVDZ ( 10) 




21- 30 


VAR I (10) 


71- 80 


D2VDR2( 10) 




31- 40 


VA R I ZP ( 10) 


81- 90 


D2VDZ2( 10) 




41- 51 


VARIZP(IO) 






CIO: 


DI MENS I 


0N( 16) 








1 


RM 


9 


ZM 




2 


R 


10 


Z 




3 


RP 


11 


ZP 




4 


DRM 


12 


DZM 




5 


DRP 


13 


DZP 




6 


DRMP 


14 


DZMP 




7 


DR PM 


15 


DZPM 




8 


DDRPM 


16 


DDZPM 


Cll : 


DIMENSI 


ON ( 56 ) 








1- 4 


BEZM ( 4 ) 


29- 32 


BEP (4) 




5- 8 


BE RM ( 4 ) 


33- 36 


BEPRP (4 ) 




9- 12 


BE ( 4 ) 


37- 40 


BEP ZP ( 4 ) 




13- 16 


BERP (4 ) 


41- 44 


DB EDR ( 4 ) 




17- 20 


BEZP ( 4 ) 


45- 46 


DB EDZ ( 4) 




21- 24 


BEPZM ( 4 ) 


49- 52 


DB EP DR ( 4 ) 




25- 28 


BEPRM ( 4) 


53- 56 


DB EPDZ (4 ) 


C13 : 


DIMENSI 


ON ( 10 ) 








1- 10 


I DDX ( 10) 






C 1 4: 


DIMENSI 


ON ( 3 ) 








1 


BZEINT 








2 


BZEMAX 








3 


BZF.F 






C15: 


DIMENSI 


ON ( 20 ) 








1- 4 


TCZDM( 4) 


13- 16 


TC DR P ( 4 ) 




5- 8 


T CDRM ( 4 ) 


17- 20 


TC DZ P ( 4 ) 




9- 12 


TCD ( 4 ) 






C20: 


D I MENS ION ( 7 ) 








1 


NDA 


5 


NAU 




2 


NRZ 


6 


NAF 




3 


NTV 


7 


NAE 




4 


NAT 






C21 : 


DIMENSI 


ON ( 150) 








1-150 


RZ ( 1 50 ) 






C22 : 


DIMENSI 


ON ( 20 ) 








1- 10 


J DDR (10) 








11- 20 


JDDZ( 10) 
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C23: 

C25 : 
C26: 
C27 : 

C30: 



DI MENS ION ( 100) 



1 




I FLAG 


59- 73 


CB SQ ( 3 , 5 ) 


2 




NC0M23 


74- 75 


DK RR DR ( 2 ) 


3- 


8 


KKRR (2,3) 


76- 77 


DK ZZDZ ( 2 ) 


9- 


14 


DDZZ (2*3) 


78- 79 


DK RZDR ( 2 ) 


15- 


24 


KKRZ (2,5) 


80- 81 


DK RZDZ ( 2 ) 


25- 


34 


CKRZ ( 2 t 5 ) 


82- 83 


DTD XP ( 2 ) 


35- 


39 


BSQT (5) 


84- 85 


DTDXM( 2) 


40- 


54 


BSQT (3,5) 


86 


DRODXP 


55- 


58 


CBSQ T ( 4) 


87 


DRODXM 








88-100 


BLANK 


DIMENSI 


ON ( 40 ) 






1- 


10 


CVRMZP(IO) 


21- 30 


CVRMZMt 10) 


11- 


20 


CVRPZP(IO) 


31- 40 


CVRPZM(IO) 


D I MENS I ON ( 16 ) 






1- 


4 


BERMZP ( 4 ) 


9- 12 


BERMZM ( 4) 


5- 


8 


BERP ZP ( 4 ) 


13- 16 


BE RP ZM ( 4 ) 


DIMENSI ON ( 32 ) 






1- 


4 


TCRMZP ( 4 ) 


17- 20 


TDRM ZP ( 4 ) 


5- 


8 


TCRP ZP ( 4 ) 


21- 24 


TD RP ZP ( 4 ) 


9- 


12 


T C RMZM ( 4 ) 


25- 28 


TDRMZM ( 4 ) 


13- 


16 


TC RP ZM ( 4 ) 


29- 32 


TDRP ZM (4 ) 


DIMENSION! 5) 






1 




INP 


4 


I NP2 


2 




I OUT 1 


5 


I0UT4 


3 




I OUT 2 
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TABLE IV 



EQUIVALENCE REFERENCE OF IMPORTANT VARIABLES 



DA ( 1 ) 




NDIM 


DA( 66) 




RM AX 


DA ( 2 ) 




IVERS 


D A ( 67 ) 




INTCON 


DAO) 




MAT 


D A ( 68 ) 




JBCZZ ( 1 ) 


DA ( 4 ) 




ITC 


DA ( 69) 




JBC ZZ ( 2 ) 


DA { 5 ) 




IDATE 


DA ( 70 ) 




BLANK 


DA ( 6 > 




JBCR( 1) 


DA ( 7 1 ) 




IDA(l) 


DA ( 7 ) 




J BC R ( 2 ) 


DA ( 72 ) 




IDA ( 1 ) 


DA ( 8) 




NDT 


# 




# 


DA ( 9 ) 




NDR 






# 


DA< 10) 




NDZ 


DA( 149) 




I D A ( 40 ) 


DA( 11) 




DT 


DA ( 1 50 ) 




I DA ( 40 ) 


DA ( 12) 




DRMIN 


DA( 151) 




BLANK 


DA( 13) 




DZMIN 


DA( 152) 




BLANK 


DA (14) 




IPIVOT 


DA ( 153) 




BLANK 


DA ( 1 5 ) 




ISTEPO 


DA( 154) 




3LANK 


D A ( 16) 




NV 


DA( 155) 




NDIMO 


DA( 17) 




I FROM 


D A ( 156) 




J J DDR 


DA( 18) 




STRTYM 


DA ( 157) 


• • 


JJDDZ 


DA( 19) 




IFYL 


DA ( 158) 


• • 


NTCDRV 


DA (20) 




NADDA 


DA ( 159) 


• • 


LASTDT 


DA ( 2 1 ) 




I IV 


DA ( 160) 




BLANK 


DA ( 22 ) 




I I VV 


DA( 161) 




BLANK 


DA ( 23 ) 




ITAPE 


DA ( 162) 




BLANK 


DA ( 24) 




IBE 


DA( 163 ) 




BLANK 


DA ( 25 ) 




RMIN 


D A ( 164) 




ICOORD 


DA ( 26 > 




ADDA(1 ) 


DA( 165) 




ISO 


DA ( 27 ) 




ADDA ( 2 ) 


DA( 166) 




JBCRR ( 1 ) 






# 


DA ( 167) 




JSCRR ( 2 ) 








DA ( 168) 




J BCR ( 1 ) 


DA (51) 




ADDA ( 26 ) 


DA( 169) 




J8CR ( 2 ) 


DA ( 52 ) 




ADDA (27 ) 


DA ( 170) 




J J DDZ ( 1) 






MATALT 


DA( 171) 




J J DDZ ( 2 ) 


DA ( 53) 




ADD A ( 28 ) 


D A ( 172) 




I SHOCK ( 1 ) 






I IAV 


DA ( 173) 




I SHOCK ( 2 ) 


DA (54) 




ADDA ( 29 ) 


DA (174) 




BLANK 






IPLAVC 


DA ( 175) 




BLANK 


DA ( 55 ) 




ADD A ( 30 ) 
IDVDX 








DA (56) 




ADDA (31) 
NREG 


DADT ( 1 ) 




TIME 


DA( 57) 




ADDA ( 32 ) 


D ADT < 2 ) 


• • 


NNDT 






11 


DADT ( 3 ) 


• • 


I S TP NO 


DA ( 58) 




ADDA (33 ) 


DADT ( 4 ) 


• • 


INEG 




NDZ2 


D ADT ( 5 ) 




MM AT 


DA ( 59 ) 




ADDA ( 34 ) 


DADT ( 6 > 




I ITC 






RZ 


DADT ( 7 ) 




JSTPNO 


DA ( 60 ) 




ADDA (35) 


DADT ( 8) 




BLANK 




NDR2 


DADT ( 9 ) 




BLANK 


DA ( 61 ) 




ADDA (36 ) 


. 




• 






Z1 


• 




• 


DA ( 66 ) 




ADDA (37 ) 










NDZ1 


DADT (25) 




BLANK 


DA (63) 




ADDA ( 38 ) 
R1 








DA ( 64 ) 




ADD A ( 39 > 
NDR1 








DA ( 65 ) 




ADDA (40 ) 
ZMAX 
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APPENDIX G 



ARRAY DESCRIPTIONS 

DA ARRAY : This array is dimensioned to 175 in C0MM0N/C3/ 

and forms a part of file RECORD. It contains control 
constants, initial values, and data labels. Its data is 
used to start, restart, and output the program. It con- 
tains, as subsets, the arrays IDA and ADDA which are 
described below. Refer to table IV of appendix F for a 
detailed breakdown of this array. 

DADT ARRAY : This array is dimensioned to 25 in COMMON/C 3/ 

and forms a part of file RECORD. It contains control 
parameters which change with every time step. Its data 
is used to keep track of time parameters in the execution 
of the program. This array, together with the VAR array, 
is buffered out to external storage on each time step which 
is saved. Refer to table IV of appendix F for a detailed 
breakdown of this array. 

VAR ARRAY : This array is dimensioned in C0MM0N/C3/ and 

forms a part of file RECORD. Its dimension must be 

determined before compilation of the program since IBM 

FORTRAN does not allow for variable dimensioning of arrays. 

This dimension is determined from: 

NVAR=NV*NDR(NV+NDZ+3) if NDR is greater than NDZ 
NVAR=NV*NDZ (NV+NDR+3 ) if NDZ is greater than NDR 

where NVAR is the dimensionality of the VAR array, NV is 
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the number of variables to be calculated, and NDR and NDZ 
are the number of mesh points in the ' R ' and *Z' directions 
respectively. After NVAR is determined, this value must 
be entered in the VAR declaration of COMMON/C 3/ in all the 
subroutines which use COMMON/C 3/ and in the NVAR assignment 
statement near the beginning of the main program. This 
must be done before compilation of the program. Refer to 
figure 1 of appendix F to find all the subroutines which 
use C0MM0N/C3/. The VAR array contains the calculated 
values of all the variables at all of the mesh points for 
a specific time step. Together with the DADT array, VAR 
is buffered out to external storage on each time step that 
is to be saved. 

RZ ARRAY : This array is dimensioned to 150 in C0MM0N/C21/ 

and forms a part of file RECORD, Actually only a portion 
of this array is used, a block equal to NDR+NDZ . The 
values in this array are used to locate a given point in 
the R-Z plane. The first NDR words give the radius, and 
the last NDZ words give the * Z* distance. Thus a mesh 
point (i,j) is located at rj=RZ(i) and z j=RZ (NDR+j ) . This 
array is buffered out to external storage after the DA 
array. 

ADDA ARRAY : This array is dimensioned to 40 in a DIMENSION 

statement and is then equivalenced to DA(26). It contains 
the values of miscellaneous additional data words. Most of 
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these values are those of initial data read into the 



initializing subroutines. This array is dimensioned to 
40 , but the number of words should be held to 26 or less 
because other variables are equivalenced to the DA array 
beginning at DA(52) (see table of appendix F). If a 
value were assigned to ADDA (27 ) by some subroutine, then 
the value of MAT ALT would also change and could cause an 
undesired effect in the execution of the program. The 
ADDA array values are assigned in subroutines MATRIX, 
TCINIT , BEXINT, MESH, and INIT4. This array together with 
the array IDA are buffered out to the line printer by 
subroutine RUNDAT. 

IDA ARRA Y: This array is dimensioned to 40 in a DIMENSION 

statement and is then equivalenced to DA(71). It contains 
the Holerith description of the values contained in the 
ADDA array. The CD C FORTRAN allows a Holerith assignment 
of eight characters to an integer variable. Thus a direct 
equivalence to the DA array could be made in that language. 
However, the IBM FORTRAN does not allow direct Holerith 
assignment to a variable. Such an assignment must be 
made through a DATA statement. In addition, the IBM 
machine will only allow four characters to be assigned to 
a REAL*4 variable. This necessitated declaring the IDA 
array to be REAL*8 and expanding the DA array by forty 
words to accommodate the equivalenced IDA array. Since 
the RUNDAT subroutine outputs the IDA array instead of 
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the DA array, no additional indexing was necessary to 
handle the two-word strings. The IDA array Holerith 
strings are assigned in subroutines MATRIX, TCINIT, BEXINT, 
MESH, and INIT4. This array together with the array ADDA 
is buffered out to the line printer by subroutine RUNDAT. 
NADDA is the total number to be printed. 

A, B, C, and V ARRAYS t The A, B, and C arrays are 10 by 10 
matrices and the V array is a ten component vector. These 
arrays are declared in COMMON/C?/ and are used in the 
difference equations 3.2-1, 3 . 5 - 1 , and 3*5-2 of UCRL 5 11 03« 

E and F ARRAYS ; These are ten component vectors which are 
declared in a DIMENSION STATEMENT and then equivalenced 
to the VAR array. They are used in forward-backward (or 
backward-f orward ) passes to solve the difference equations 
in sections 3.2 and 3*5 of UCRL 5H03. 

W ARRAY ; The W array is a reference vector which con- 
tains all the variables, excluding alternate variables, 
which can be calculated by this program. Alternate 
variables are derived from selected W vector components 
by multiplying or dividing by the density. The vector 
W is defined by W=(R 0 ,VR,VP ,VZ,TI ,TE,BR,BP,BZ ,SI) . 

U ARRAY : The U array is a rearrangement of the W vector. 

The order of the rearrangement is determined by the IV, 

IIV, and IAV transformation vectors which in turn are 
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determined by JAV and JVERS in subroutine MATRIX . After 
U is determined, only the first NV components are calcu- 
lated, NV also being automatically determined in sub- 
routine MATRIX. JAV and JVERS are parameters on input 
card number 7. 

TEMP ARRAY : The TEMP array is a temporary vector identical 

to the U vector. It is used during calculations and con- 
tains the values of the U vector along a single line 
adjacent to the U vector line at the new time step. For 
example, during calculations the U array contains U n+1 for 
line *k'. At that time, the TEMP array would contain U n+1 
for line • k-1 • . 

IV, IVV , and IAV ARRAYS : These arrays are determined 

automatically by JAV and JVERS in subroutine MATRIX. The 
arrays are assigned values from integer constants IIV, 

II W, and I IAV which are assigned in the subroutine. 

IV(1) is the value of the first digit of IIV, IV(2) the 
value of the second, and so forth* The IV array determines 
the arrangement of the U vector. This is best explained 
by an example. If IV (6) =4, then the sixth component of 
W (TE) is placed as the fourth component of U. Again 
if IV(3)=6, then the third component of W becomes the sixth 
component of U, NV determines the number of components 
which will actually be calculated. Only the first NV 
components of U are calculated. Therefore, if IV (K) is 
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greater than NV, then the Kth component will not be 
calculated. 

The IVV array essentially gives the same information 
that the IV array does, only in reverse. For example, 
if IV(6)=4, then IIV(4) would equal six and would in- 
dicate that the fourth component of U is the sixth component 
of the reference vector W. 

The IAV array operates in a similar manner and is 
concerned with alternate variables. Here IAV(6)=1 says 
that the first alternate variable for the sixth component 
of W is to be calculated. In this case FE (electron 
pressure) is calculated instead of TE (electron temperature). 
One further example should put these transformations in 
perspective . 

W=(RO ,VR ,VP ,VZ ,TI , TE ,BR ,BP , BZ ,SI ) 

I IV= 1267348950 with NV=5 

IIVV= 125693^780 
U=(R0. VR. TI. TE.BZ. VP.VZ.BR, BP t Si ) 

’not calculated since NV=5 

These arrays are declared and passed to other subroutines 
in C0W0N/C6/. 

Use of these transformation arrays permits efficient 
utilization of storage and time. The code is written to 
calculate all of the variables in W. If not all of these 
variables are required, much time and storage area would 
be wasted in their calculation. These arrays allow the 
operator to select the variables to be calculated and 
suppress all the others. Several subsets of W are already 
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set up in subroutine MATRIX. Any additional subsets which 
may be required would have to be added to MATRIX. The 
procedure for adding subsets is outlined in the description 
of subroutine MATRIX in appendix D. 

IBCRR and IBCR ARRAY S; These arrays are dimensioned 
and declared in COMMON/C 4/ and determine the inner and 
outer *R' boundary conditions respectively. The array 
values are computed from the values of integer constants 
JBCRR and JBCR which are read from input card number 5A. 
IBCRR(l) is the value of the first digit of JBCRR, IBCRR(2) 
is the value of the second digit, and so forth. The same 
procedure determines the values for the IBCR array. The 
following example demonstrates the use of these arrays 
in determining the boundary conditions. If JBCRR=3127463100 , 
then IBCRR ( 1 )=3 • IBCRR(2)=1, IBCRR(5)=4, and so forth. 

IBCRR( 1 )=3 indicates that the first component of W (RO) 
has inner radial boundary condition number three. Again 
IBCRR(5)=4 indicated that the fifth component of W (TI) 
has inner radial boundary condition number four. Boundary 
conditions are discussed under the subroutine BCR in 
appendix D, 

IBCZZ and IBCZ ARRAYS : These arrays are dimensioned and 

declared in COMMON/C 4/ and determine the lower and upper 
' Z* boundary conditions respectively. The array values 
are computed from the values of integer constants JBCZZ 
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and JBCZ which are read from input card number 5A. The 
array values are determined and utilized in the same manner 
as the IBCRR and IBCR array values are used. Boundary 
conditions are discussed under the subroutine BCZ descrip- 
tion in appendix D. 

JDDR and JDDZ ARRAYS : These arrays are declared and 

dimensioned in C0MM0N/C22/ and determine the use of one- 
sided derivatives (page 79 of UCRL 5*103 ). The array 
values are computed from the values of integer constants 
JJDDR and JJDDZ which are read from input card number 6. 
This computation scheme is the same one used for IBCRR 
and IBCR above. The array elements can only take on the 
values '1* or 'O’. Zero indicates no one-sided derivative 
and one indicates use of the one-sided derivative. If 
JJDDR(5)=1» the fifth component of W (Tl) will be cal- 
culated using one-sided derivatives. 

NOTE: Throughout this paper the terms array and vector 

are sometimes used interchangably . This arises from the 
changes Lindemuth made to conserve storage area. The 
variable arrays were rewritten as singly subscripted vari- 
ables and thus became vectors instead of arrays. In the 
context of this paper the vector may be considered to be 
a one -dimensional array. 
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APPENDIX H 



DATA SETS 

A data set is a collection of information which 
resides on a particular device. Typical data sets are 
the card reader, card punch, line printer, plotter, 
magnetic disc, magnetic tape, data cells, and paper tape. 

The magnetic core in the CPU is also a data set. It is 
an internal data set and is not treated like the external 
data sets mentioned above. Data sets must be declared 
and defined in the JCL program. 

The data set is declared and defined in the data 
definition (DD) statement. The DD statement is composed 
of three sections which are (from left to right) the ident- 
ification field, the operation field, and the operand field. 

A. IDENTIFICATION FIELD 

In this field the JCL program stepname and data set 
identifier are specified. The following identif ication 
field 

//G0.FT10F001 

identifies a data set FT10F001 which is used in the GO 
step of a program or catalogued procedure. The third and 
fourth digits of this identifier form the integer device 
identifier which is used in the READ/WRITE statements in 
the source program. For example READ (10 ,50) A would cause 
variable 'A' to be transferred from core to the device 
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which contains data set FT10F001. 

The computer facility at the Naval Postgraduate School 
has specified that data sets FT05F001 , FT06F001, and 
FT07F001 reside on the card reader, line printer, and 
card punch respectively. These data sets do not have to 
be defined in the JCL program. If these devices were to 
be used with different data set identifiers, defining DD 
statements would have to be included in the JCL program. 

Data sets FT06F001 and FT07F001 are used in the GO 
step of the common catalogued procedures. If the card 
punch or line printer is to be used during the FORT step 
(compilation of FORTRAN) the data sets would be SYSPRINT 
and SYSPUNCH respectively. These data sets do not have to 
be declared by the operator. Frequently SPACE parameters 
in these data sets must be over-ridden when a very large 
listing or object deck is expected during compilation of 
a source deck. The SPACE parameter is discussed in para- 
graph C below and over-ride procedures in appendix B. 

B. OPERATOR FIELD 

The operator field in all data definition statements 
contains 'DD'. This identifier must be preceded and followed 
by at least one space. 

C. OPERAND FIELD 

This field contains all the information necessary to 
completely describe the data set. The following parameters 
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follow the DD entry and may appear in any order. However, 
they must be separated by a comma and no space. Spaces 
are not permitted to appear anywhere within the operand 
field except after the end of the field or on a continuation 
card between the comma following the last parameter and 
the ' C 1 in column 72. A space anywhere else will termi- 
nate the program. 

1 . SPACE Parameter 

This parameter specifies the amount of storage 
to be allocated for the data set. Space is allocated by 
TRACKS (TRK ) or CYLINDERS (CYL). A track can store 7294 
bytes of information. A cylinder consists of 20 tracks. 

For general information, a 2314 disc contains 200 cylinders. 

The space parameter size is determined from the 
data storage requirement and the BLOCKSIZS specified in 
the DCB parameter discussed below. First determine the 
total amount of storage in bytes required to hold the expected 
output data. Divide this figure by the blocksize to de- 
termine the total number of blocks in the data set. Next 
divide this figure by the largest number of blocks which 
will fit in a track. This will give the number of tracks 
required. If allocation is to be made in cylinders, divide 
the number of tracks by 20 and add one for any fractional 
part of a cylinder. 

Allocation in tracks allows the machine more flex- 
ability in locating space since the tracks do not have to 
be consecutive. However, there is advantage to allocating 
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space in cylinders. If it is anticipated that a large 
portion of the data set will be used at one time, alloca- 
tion should be in cylinders because then the mechanical 
read/write arm does not have to move during operation on 
a particular cylinder. If the allocation were in tracks, 
it is conceivable that the arm would have to move after 
the operation on each track was completed. This could 
be very time consuming if many i/o commands were made to 
that data set. 

The SPACE assignment will take the form 
SPACE* (TRK, (15.2) .RLSE) 

which will allocate 15 tracks initially. If the 15 tracks 
are insufficient, additional tracks will be allocated two 
at a time until 18 additional tracks have been added. RLSE 
will release any unused area after the program has termi- 
nated. The RLSE parameter is normally used only if the 
required size of the data set is unknown and a very high 
estimate is used to assure enough space. Thus the SPACE 
assignment frequently appears in the form 
SPACE=( CYL , (3.1)) 

where three cylinders are initially allocated with incre- 
ments of one being allocated as needed. Care must be 
taken that all commas and parenthesis appear in the proper 
location. This has been a source of numerous errors on the 
part of the author during this conversion. 
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2 . UNIT Parameter 



The UNIT parameter specifies the device upon 
which the data set will reside. There are identifiers 
for each device in use with the CPU. This paper will be 
concerned only with magnetic discs for data sets other than 
the card reader, line printer, and card punch. 

For permanent data sets, the UNIT parameter will 
always be 2314 which is the identifier for a magnetic 
disc device which is used by general users. There are 
presently three discs available for the general user. The 
particular disc on which the data set is to reside is speci- 
fied by V 0L- SE R= * name 1 where 'name* is either MARY, LINDA, 
or DUFFY. The operator should check with the dispatcher 
to assure that adequate space is available for a particular 
disc before completing this parameter. The UNIT assignment 
will then take the form 
UNIT=2 314, VOL=SSR=MARY 

which places the data set on a 2314 disc pack named MARY. 

For temporary data sets, the assignment is UNIT= 
SYSDA. SYSDA is reserved for general user temporary data 
sets. The data sets are destroyed after the program ex- 
ecution is terminated. When UNIT=SYSDA is used, VOL=SER= 
'name' and DISP, which is discussed below, are not required. 
This unit assignment was used for all of the data sets during 
the testing phase of the conversion when it was not im- 
portant that the information be saved. 
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3 . Data Control Block (DCB) Parameter 

The DCB parameter describes the construction of 
the data set. It is composed of three parts; the record 
format, the logical record size, and the block size. 

a. Record Format 

The record format (RECFM) specifies the type 
of grouping of data elements in the data set. Typical 
RECFM' s are fixed blocked, variable blocked, variable 
blocked spanned, and so forth. This program uses fixed 
blocking and the assignment if RECFM=FB. 

b. Logical Record Length (LRECL) 

The logical record size is the length in 
bytes of the list in the READ/V/RITE statement used for 
the data set. For example, the statement 
WRITE (10) A , B ,C 

where A, B, and C are REAL*4 variables would require a 
logical record length of 12 since the list contains three 
words at four bytes per word. There are eight bits per byte. 

In this program, file RECORD is placed in data 
set FT10F001. The LRECL for this data set is 100. Thus 
to read or write vectors of greater length that 25 words, 
the i/o commands were rewritten in DO loops so that each 
READ/WRITE statement has a list of only 25 words in length. 

c. Block Size (BLKSIZE ) Parameter 

The block size parameter is specified to be 
an integer multiple of the LRECL. The normal practice is 



96 



to dimension the BLKSIZE such that two blocks will fit on 
one track. Thus the maximum BLKSIZE will be 3520. To 
determine the block size for the data set, divide 3520 
by the LRECL. The truncated integer result multiplied by 
the LRECL will be the BLKSIZE in bytes. In this program 
the BLKSIZE is 3500. 

The DCB parameter will take the form 
DCB=(RECFM=FB , LRECL=100 , BLKSI ZE=3500 ) 
which indicates that the data set consists of fixed blocks 
of 35 logical records. The logical records are 100 bytes 
long and there are two blocks per track of disc. 

4. Disposition (DISP) Parameter 

The disposition parameter defines the status and 
disposition of the data set and is composed of two sub- 
fields. 

a. Status 

This sub-field may contain either OLD, NEW, 

MOD, or SHR. OLD indicates an existing data set while NEW 
indicates a data set that is being created. MOD indicates 
an existing data set which may be modified and SHR indi- 
cates that the data set may be shared by other jobs. 

b. Disposition 

This parameter may be either DELETE, KEEP, or 
PASS. DEIETE means that the storage occupied by the data 
set will be released, when the job is terminated. KEEP says 
to keep the data set for later use, and PASS says to pass 
the data set to succeeding steps in the job. 
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Thus the statement DISP=( NEW .KEEP) means keep 
the data set being created after the job is over. 

5. Data Set Name (DSNAME) Parameter 

DSNAME names the data set and takes the form 
DSNAME=Unnnn.anyname where •U* indicates the type of user, 
nnnn is the user’s assigned number, and ’anyname ' is a unique 
name assigned by the user to identify the data set. This 
name is limited to six alphanumeric characters. The type 
of user will be 'S' for student, *F* for faculty, 'C' 
for computer center staff, * N • for NFS staff, and 'X' for 
external users. The statement DSNAME=S0833 . RECORD gives 
the name RECORD to a data set used by a student v/hose 
user number is 0833. The DSNAME parameter is not required 
when defining temporary data sets on unit SYSDA . 

6 . Label (LABEL) Parameter 

The label parameter is used to indicate the ex- 
piration date of a permanent data set and may take two 
forms : 



LABEL=EXPDT=yyddd where yy is the year and ddd is the 
Julian day of the year. 

LABEL=RETDP=zzz where zzz is the number of days from 
its creation that the data set is to be retained. 

The statement LABEL=EXPDT=?2366 indicates that the data 

set will be retained until 3^ Dec 72, 



A complete DD description will now take the form 

//GO . FT10F001 DD UNIT=23l4 ,DSNAME=S0833 . RECORD 
// DIS?=( NEW, KEEP) , V 0L= SE R= M A R Y , 

// LABEL=SXPDT=?2366 t SFACE=(TRK , (20.2) ) , 

DCE= ( RECFK=FB , LRE:L=100 , BLKSI ZE=3500 ) 



2 



c 

c 

c 
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The source program statement WRITE (20) A will transfer 
the 32 bit pattern under the identifier 'A* from core into 
the data set named RECORD. This data set is located on 
disc pack MARY and will be retained until 31 Dec 72. 
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APPENDIX I 



IDENTIFIER DESCRIPTIONS 

This section contains brief descriptions of the initial 
value and control identifiers. The following is not a 
complete list of all the identifiers used in the program, 
but rather a list of the identifiers which set initial 
values and control the execution of the program through 



subroutine 


argument lists and COMMON blocks. 


BP 


Component of B Held in the phi direction 


BPF 


Frequency plasma B field in the phi direction 


BPINIT 


Initial plasma B field in the phi direction 


BPMAX 


Maximum plasma B field in the phi direction 


BR 


Component of B field in the R direction 


BRAT 


Mirror ratio 


BZ 


Component of B field in the Z direction 


BZF 


Frequency of plasma B field in the Z direction 


BZEF 


Frequency of external B field in the Z direction 


BZEINT 


Initial external B field in the Z direction 


BZEMAX 


Maximum external B field in the Z direction 


BZINIT 


Initial plasma B field in the phi direction 


BZMAX 


Maximum plasma B field in the Z direction 


CLNLM2 


Constant used in modified Coulombian logarithm 


CQA 


Shock smoothing used 


DT 


Size of time step 


ECHO 


Electron charge 



100 



EMASS 


Electron mass 


EO 


Permittivity 


GAM 


Gamma Factor 


I ARRAY 


Selects array to be buffered in or out by BUFFER 


IBE 


IBE=1 separates B field into plasma and external 
components and reads in external field components 


IDATEO 


Julian data, five character 


IFROM 


Selects restart mode. IFROM/O restarts from 
file OLDREC. IFRCM=0 restarts from file RECORD. 


IFYL 


File reference number 


INP 


Data set FT05F001 identifier 


INP2 


Data set FT30F001 identifier 


INTCON 


Selects initial condition subroutines - 
* 1 • selects INIT1 , ’4’ selects INIT4, and 
•5* selects INIT5. 


I OP 


Selects profile option 


IOUT1 


Data set FT10F001 identifier 


I0UT2 


Data set FTOoFOOl identifier 


I0UT4 


Data set FT40F001 identifier 


IPU 


Selects restart mode - *0' allows new parameters 
to be read in MATRIX on restart. *1' causes 
code to use original parameters 


IROR 


Determines read or write function in BUFFER 


ISHOCK 


Selects variables to receive shock smoothing 


ISO 


Selects isotropic or anisotropic thermal 
conductivities 


1ST 


Gives the number of the time step from which 
restart is to begin 


ISTEP 


Selects coding for different passes through 
equations 


ISTEPO 


Used to set ISTE? odd or even 
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ITAPE 


Tape reference number 


ITCVER 


Subroutine version 


ITP 


File identifier in READ/WRITE and BUFFER commands 


JAV 


Selects alternate variables to be calculated 


JBCR 


Sets outer R boundary conditions 


JBCRR 


Sets inner R boundary conditions 


JBCZ 


Sets upper Z boundary conditions 


JBCZZ 


Sets lower Z boundary conditions 


JDDR 


Selects one-sided derivatives for R 


JDDZ 


Selects one-sided derivatives for Z 


JVERS 


Selects variables to be calculated 


J2BCHK 


Gives R coordinate of mesh point to be checked 


K2BCHK 


Gives Z coordinate of mesh point to be checked 


MAT ALT 


Selects coding subset in MAT2, This code uses 
MATALT=0 . 


MMAT 


Selects difference method 


NDIM 


Number of dimensions 


NDR 


Number of mesh points in the R direction 


NDT 


Number of time steps 


NDTBUF 


Arrays DADT and VAR will be saved every NDTBUF 
time step 


NDTPNT 


0UTPT1 will be called every NDTPNT time steps 


NDZ 


Number of mesh points in the Z direction 


NREG 


If greater than ' 1 ' , card 11B is read 


NTCDRV 


NTCDRV=1 indicates no transport coefficient 
derivatives 


NZP(K) 


Number of the Kth Z line to be printed by 0UTPT1 
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N2BCHK 


Gives number of mesh points which will be 
checked by CHECK and SPECHK 


PMASS 


Mass of the proton or ion 



PI 



RMAX 


Maximum radius 


RMIN 


Minimum radius 


RO 


Density 


ROINIT 


Peak density 


ROMIN 


minimum density 


RRAT 


Determines mesh spacing factor in the R direction 


SI 


PSI, the stream function 


SMASS 


EMASS+PMASS 


STKEQO 


STKEQC= 1 sets thermal conductivity equal to zero 


SV 


Specific volume 


TE 


Electron temperature 


TI 


Ion temperature 


TINIT 


Peak temperature 


TMIN 


Plasma background temperature 


UO 


Permeability 


VEQMAX 


Where velocity is maximum 


VEQO 


Multiple of XO where velocity goes to zero 


VO 


Plasma initial boundary velocity 


VP 


Component of velocity in the phi direction 


VR 


Component of velocity in the R direction 


VZ 

XO 

XO 


Component of velocity in the Z direction 
Plasma initial boundary radius 
Plasma initial boundary radius 
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ZMAX 



Maximum distance in the Z direction 



ZRAT 



Determines the mesh spacing factor in the 
Z direction 
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CONVERTED SOURCE PROGRAM 3 JUN 72 
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